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ABSTRACT 


Major results of a two year project to study the climatological and micro- 
meteorological conditions in a tropical evergreen forest region are summarized 
in this report. The study has been under the sponsorship of the U. S. Army 
Engineer Topographic Laboratory, The forest area under study is in the 
interior of Thailand and is influenced by a cool, dry northeast monsoonal 
flow from November to March and a warm, moist southwest monsoonal flow 
from May to September. Measured information were collected during a three 
year field program sponsored by the U. S. Army Natick Laboratories and 
included ground station temperature, relative humidity, precipitation, 
evaporation, hours of sunshine, and tower measurements of temperature, dew 
point temperature, and wind speed, incoming and reflected solar radiation, 
incoming and outgoing infra-red radiation and sub-surface temperature 
profiles in the forest and in a cleared area within the forest region. ^ 

The results discussed here include a climatological survey of the 
experimental site as well as a comparison of the monsoonal cycles at a 
number of widely dispersed sites in Southeast Asia; an analysis of the 
profiles of temperature and wind speed in the forest and cleared area; 
diurnal variations of static stability and Richardson number representative 
of several different periods of the year; a summary of calculations of the 
roughness characteristics of the forest canopy and the variability thereof; 
spectral analyses of kinetic energy above and within the forest canopy 
during both monsoonal flow regimes; estimates of the albedo of the forest 
and the net fluxes of radiational energy for several periods of the year; 
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SECTION I - INTRODUCTION 


This Technical Note serves as a final report on Contract DAAK-02-72- 
C-0287 from the U. S. Army Engineer Topographic Laboratory, Ft. Belvoir, 
Virginia, to the Graduate Program in Meteorology, University of Maryland. 

The scope of the research under this contract was to prepare a scientific 
analysis of the climatological and micro-meteorological data obtained during 
a scientific field program called TREND (Tropical Environmental Data). 

The TREND field experiment was under the sponsorship and management of 
the Earth Sciences Laboratory, U. S. Army Natick Laboratories (Dr. Paul 
Dalrymple: Project Supervisor) and carried out by the Applied Scientific 
Research Corporation of Thailand in collaboration with several agencies 
of the Thai Government. The experiments were conducted during the years 
1967-70 as an interdisciplinary study of a tropical dry evergreen forest 
environment and embodied surveys of micrometeorological conditions; 
hydrologic factors; soil, vegetation, and biological surveys. Meteorological 
data from the TREND experiments, which serve as the source of information 
for the scientific analysis in this report, were provided by USAETL in tlie 
form of raw data tapes recorded at the TREND site and certain computer 
printouts prepared during the field experiment. 

Personnel participating in the processing, analysis and interpretation 
of the TREND meteorological data were as follows; Dr. Helmut E. Landsberg, 
Principal Investigator; Dr. Owen E. Thompson, Co-Investigator; Mr. Robert 
E. Kay lor. Computer Systems Specialist; Mrs. Rachel T. Pinker, Mr. David 
K. Lee, Dr. Ilassan El-Sayed, Miss Jae K. Eom, Project Research Assistants; 

Mr. Mark Rodney, Technician; Mrs. Syvia Epstein, clerical assistant. 
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The study site for the TREND experiment is the Sakaerat Forest, an area 
of about 80 square kilometers in Thailand situated on the Korat Plateau at 
approximately 1A°31'N, 101°55'E which is about 60 km south-southwest of 
Khorat ( Nakhon Ratchasima) and about 190 km northeast of Bangkok. (See 
Figure 1) Meteorological data were taken along two towers approximately 
50 ra high and 450 m apart, one tower penetrating the forest, the other 
erected in a small clearing and both based at about 535 m above sea level 
on a tilted slightly dissected sandstone plateau which slopes from southwest 
to northeast* Figure 2 shows the forest area under consideration and the two 
towers mentioned above. The forest Includes good stands of dry evergreen 
forest, dry dipterocarp forest and clearings in various stages of regenera- 
tion. The immediate vicinity of the towers Us comprised oj a two-storied 
forest. The top story consists almost exclusively of llopea ferrea , with a 
scattering of Shorea sericeif lora . with dense and continuous crown canopy 
with tree tops ranging from 20-35 m high. The second story consists of the 
species Hydrocarpus iliclfolius . Walsura trichostemon , Aglaia pirlfera . 
Lagerstroemia . and Memecylon ovatum with tops ranging from 5-17 m high. 

The forest floor was well covered with thick undergrowth. 

The experimental design for the micro-meteorological survey called 
for the collection and processing of temperature, dew point temperature, 
wind speed and direction at numerous levels along both towers, various 
radiation fluxes at the tops of the towers and near tlie ground, and sub- 
surface temperature profiles near both towers. Six satellite ground stations 












recorded temperature, humidity, precipitation and potential evaporation. 

Two base stations, each within about 2.5 km of the towers, also had sunshine 
recorders, additional rain gauges and evaporation gauges. Figure 3 shows, 
schematically, the location of the various measuring Instruments along 
each tower. A general account of the instrumentation and monitoring 
procedures is given in ASRCT (1969). Briefly, the data logging system 

ff y ‘X » 

[Hewlett Packard] e omp? ioc d - two'systems combined into one, capable 

A 

of analogue and digital recording. The digital system records the 42 values 
of temperature and 16 values of wind speed, along with date and time. The 
analogue system records the various radiation measurements, wind direction, 
precipitation, dew point temperature and station pressure. The frequency 
of scanning in both systems is of the order of once per half-minute. In 
what is to follow, data recorded on the digital system is referred to as 
"D-tape data", that on the analogue system as "A-tape data". 

The harshness of the tropical climate — heat , moisture , mildew, etc. — 
made maintenance of the measuring systems very difficult. In particular, 
the dew point sensors frequently malfunctioned so that the accuracy and 
reliability of dew point temperatures in the data record must always be 
suspect . 

The results presented in this report emanate from efforts at the 
University of Maryland to extract useful climatological and micro-meteoro- 
logical information from the raw data available to us from Project TREND. 
Section II serves as an outline of the extensive development of data reduc- 
tion schemes which were necessary to convert the raw data material into a 


temperature: all levels 

DEW POINT: ALL LEVELS EXCEPT 16 

WIND SPEED, direction: LEVELS WITH SOLID LINES 

precipitation: levels 1 . 16 

SUBSURFACE TEMPERATURE: 5 . 10 . 20 . 50.100 CM 

radiation: 

FOREST CLEARING 

LEVEL 1 : RjCdIRECT) LEVEL 14 : R^^ -- Rg'^ 

LEVELS 1 . 14 : R|_* . R^-^ 

Rg* . Rg-f . R^(SPECTRAL) 

Ri(UV) . ILLUMINATION 




FOREST CLEARING 
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data resource useful for analysis. Section III sununarizes the climatological 
conditions at the TREND experimental site while at the same tine placing 
this into perspective with regard to the monsoonal cycle in Southeast 
Asia. Section IV summarizes results of a careful analysis of the wind 
speed and temperature distributions in the forest environment. Section V 
summarizes the behavior of various energy fluxes in and above the forest 
canopy. 
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SECTION II - PROCESSING OF TREND EXPERIMENT 
RAW DATA AI'ID CURRENT DATA AVAILABILITY 

A considerable fraction of the effort during this project was devoted 
to tlie preparation of the raw micro-meteorological data for computer 
processing. This proved to be a formidable task because of the condition 
of the raw data tapes and the manner in which data are recorded on tliose 
tapes. Computer techniques were developed for the reading of raw "D" 
tapes, on which the basic wind and temperature profile data are recorded, 
and repacking the information on magnetic tape in a form convenient for 
liigh speed processing on a Univac 1108 machine but also compatible with 
Control Data Corporation hardware. The following discussion summarizes 
the procedure. In the discussion, the term "word" refers to a 36 bit, 
Univac 1108 computer word. 

On the raw D-tapes from the TREND Field Experiments, the time, wind 
speed and temperature data were recorded every 10 seconds in a 'BCD' format 
on an incremental tape recorder. Each scan consists of 1 digit for the 
year, 3 digits for the day and 2 digits each for the hours, minutes and 
seconds. Following this are A2 temperature points each consisting of a 
sign, channel number, and 4 digits of temperature. The last group consists 
of 3 digits of wind speed, every other one starting with the letter ' F' . 

One scan requires 102 words or 612 6-bit bytes. There are 1024 words per 
tape record, (6144 bytes), or roughly 10.04 scans per record. At 556 
bytes per inch recording density, plus 0.75 inches for the inter-record 
gap, each record is about 11.8 inches long. The major problem in reading 




these Capes is that the incremental recorder used tended to stretch the 
tape, resulting in poor spacing between characters. On a modern tape 
1 drive system (Univac, IBM, CDC, etc.) this spacing is interpreted as an 

inter-record gap, resulting in loss of data or loss of position on the tape 
drive. The amount the tape was stretched can be seen from the fact that 
one of the unreadable tapes has 2739 records, which requires almost 2700 
feet of tape, but a normal reel only contains 2400 feet of tape! 

After considerable experimentation, it was concluded that an older 
tape drive could read these tapes better than a more advanced tape drive. 

1 The University of Maryland still owns two IBM 1401 systems which seem to 

j be well-suited for this copy work. Because the 1401 has only 8000 bytes of 

! 

j memory, and the input record requires 6144 bytes, the copy program was made 

as simple as possible, 

Tlie program functions by first clearing the core, then reading in a 

I record. The record is accepted regardless of parity errors, abnormal 

( 

frame counts, input-output checks, etc. Then 1024 words (36 bit, Univac 
1108 word) are written on a temporary tape in even parity, 800 bytes to 
I the inch. If the input record is shorter than this it is padded with BCD 

! 'Z's. If for some reason it is longer, it is chopped to this length. The 

program v'ill generally recover and recopy 90% or so of tlae data on the ori- 
ginal raw "D" tapes obtained from the TREND site recorders. 

The next step in the procedure is to repack the data from the temporary 
tapes, obtained from the 1401 processing routine, onto permanent storage 
tapes in a compressed format compatible with Univac 1108 and CDC 6600, 7600 



■j 


j 


j 

! 



viii 




hardware. These final storage tapes then represent the data base for all 
subsequent computer processing for TREND micro-meteorological studies. Tlie 
compression factor is about 5 to 1. The following is a summary of the 
program labeled 'PACK' vjliich produces permanent storage tapes from the 
temporary copy tapes. 
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PROGR^XM 'PACK' FOR CONVERSION OF TEMPER.\TURE Al.'I) WIND DATA TAPES, 

INTRODUCTION 

THIS PROGR.\M ACCEPTS FIVE RAW DATA TRAPES AS INPUT AND PACKS THEM 
ONTO ONE OUTPUT TAPE. THE RAW INPUT T/M’ES ARE IN EVEN PARITY, BCD CODE, 

.VND 556 BPI DENSITY; THE OUT-PUT TAl’ES ARE ODD PARITY, BINARY, /VND 300 BPI 
DEN’S ITY. ,UL TRAPES ARE 7-TRACK. THE PACKING CONVERTS EACH DATA ITEM TO AN 
EIGHTEEN BIT QUANTI'TY AND STORES TWO ITEMS PER COMPUTER WORD (36-BIT, 

1108 COMPUTER). 

THIS PROGRAM HAS BEEN DESIGNED SPECIALLY FOR THE UNIVAC 1108 COMPUTER; 

IN PARTICULAR, THE 18 BIT ITEM LENGTH WAS SELi;CTED TO T/\KE /ADVANTAGE OF THE 
PART-WORD TR.ANSFERS AVAILABLE ON THIS MACHINE. HOWEVER, THE FINAL RECORD 
LENGTH WAS AD.IUSTED BY PADDING SO THAT EACH RECORD IS AN EXACT MULTIPLE 
OF SIXTY BITS AND THEREFORE IS COMPATIBLE WITH THE CDC 6600 OR 7600 SERIES 
COMPUTERS . 

AN ATTEMPT HAS BEEN MADE TO RECOVER AS M/ANY OF THE B/AD INPUT RECORDS 
AS SEEM ECONOMICALLY FEASIBLE; A SAMPLE RUN WITH A POOR INPUT TAPE RECOVERED 
ABOUT 94% OF THE DATA. WHENEVER A TIME SC/AN CANNOT BE RESOLVED, THE SC.-VN 
IS 'FLAGGED' AND 'THE PROPER TIME UNITS INSERTED WITH ZERO FILL FOR IHE DATA 
UNITS. IF THE TIME G/AP IS OVER 60 SECONDS, NO FLAGGING OR FI LUNG IS DONE; 

THE RECORD CONTINUES Wmi THE NEW TIME. 

ALL ITEMS ARE CARRIED AS FIXED POINT INTEGERS. THE TEMPERATURE DATA 
HAS OCTAE 300,000 ADDED TO AELOW FOR THE POSSIBIEITY OF NEGATIVE TEMPERATURES; 
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UTIEN THIS QUANTITY IS SUBTRACTED, THE ORIGINAL DATA IS RECOVERED EX;\CTLY. 

THE WIND SPEEDS ,\RE IN PULSES PER SECOND, AND MUST BE DIVIDED BY 31.87 AND A 
CONSTANT OF 0.52 ADDED TO THE RESULT TO OBTAIN MILES PER HOUR. THE FACTORS 
3187 AND 52 ARE PACKED WITH THE DATA. NO CHANNEL NUMBERS ARE PACKED ONTO 
THE OUTPUT, BUT THE DATA IS ALWAYS l^rRITTEN IN THE SAME ORDER AS THE INPUT 
DATA: TEMPER.\TURE CILYNNEL NOS. 101, 201, 301, .... 114, 214, 314; THEN WIND 
CHANNELS 1, 2, . . ., 16. 

INPUT TAPES 

FIVE TAPES .ARE WRITTEN SERI/ALLY ONTO ONE OUTPUT T^APE. THE PROGR/VM 
ASSUMES TH.AT THESE FIVE TAPES ARE IN ORDER AND FOLLOW EACH OTHER IN TI.’^E. 

IF THE LAST RECORD ON THE FIFTH TAPE DOES NOT END WITH A COMPLETE SCAN, 

TH"S I.AST SCAN WILL PROBABLY BE LOST. THE INPUT TAPES MUST BE ASSIGNED WITH 
FORTR.AN NUMBERS 10, 15, 20, 25, 30 IN TilAT ORDER (EXEC. II UNIT NOS. "D", "I", 
"N", "S", "a"). if more than one OUTPUT TAPE IS TO BE RHITTEN, EACH I.NPUT UNIT 
MAY BE RELOADED .AS SOON AS THE PROGR^VM FINISHES WITH IT. THE PROGRAM WIIJ. 
REWIND THE TAPE UNITS AS THEY .■ARE FINISHED. NOTE IIIAT THE PROGRAM EXPECTS 
FIVE INPUT TAPE I'NITS TO PE ASSIGNED TO IT, RATHER THAN ONE MULTI-REEL UNIT. 

EACH INPUT T.APE UNIT MUST BE ASSIGNED WITH EVEN PARITY, 800 BPI DENSITY 
(FOR 1401 COPY TAPES), HARDWARE TR.ANSL/ATE TO FIELD-DATA, /AND PARITY ERROR 
OPTION IF AVAILABLE. ESTIMATED TIME PER INPUT TAPE IS BETWEEN TWO AND 
THREE MINUTES, DEPENDING ON THE NUMBER OF RECORD ERRORS. INPUT IS COMIM.ETELY 


BUFFERED SO THAT VERY LITTLE TIM! IS LOST ON T/APE READS. 
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t)UTPUT 

OUTPUT TAPES ARE ASSKINED FORTRAN NUMBERS 12 ANU 14, IN THAT ORDER. 

(EXEC. II UNITS "F" AND "H") . IF Mt)RE THAN TWO OUTPUT TAPES ARE TO BE 
WRITTEN, THE TWO UNITS ARE RELOADED AS NECESSARY. THE PROCIUM WILL REWIND 
EACH UNIT AS IT IS FINISHED. 

AGAIN, THE PROGILXM EXPECTS TWO PHYSICAL UNITS TO BE ASSIGNED, RvYTHER 
THAN ONE MULTI-REEL UNIT. THESE UNITS MUST BE ASSKINED WITH ODD PARI'n’ , 

800 BPI DENSITY, NO TRANSLATE. EACH OUTPUT INYPE REQUIRES ONE CARD AS DATA 
INPUT; THIS IS THE REEL NUMBER AND MUST BE A THREE DIGIT INTEGER BETWEEN '001' 
AND '511' (OCTAL 777). (FORTR/VN FORMAT 13). THERE MUST BE ONE DATA CARD 
PER OUTPUT TAPE FOLLOWED BY ONE CARD WITH '000' TO STOP THE PROGRAM. IF 
THIS LAST CARD IS LEFT OUT, THE PROGRAM WILL EXIT BY 'ERRS' RATHER THAN THE 
NORMAL EXIT, BUT OTHERWISE NOTHING DR^XSTIC WILL HAPPEN. 

LOAHNG 

'PACK' IS AN ASSEMBLY LANGUAGE PROGR.\M WHICH SHOULD COMPILE WITH ANY 
1108 ASSEMBI.Y AINGUAGE COMI’ILER. IT REQUIRES LIBRy\RY ROUTINES 'NWDUS', 
'NTRAN', 'NSTOP$', 'NRDU$', 'NIOl$', 'NI02$', AID 'NINTHS'. IF THIS LAST 
ROUTINE IS NOT AVAILABLE, IT M/\Y BE ELIMINATED BY DELETING CAIDS 98 AND 
99 WrniOUT AFFECTING THE PIUXIR/AM. IF ANY OF THE OTHER ROUTINES ARE l.,VBELED 
DIFFERENTLY, REFER TO THE REFERENCE TABLE AND CHANGE THE PROPER CARDS. 

THE PR0(;R.\M also RE()UIRES one PROCEDURE, 'AXR$' WHICH DEFINES REGISTERS, 
IMMEDIATE OPERATIONS, PART WORD TR/\NSFERS, ETC. THIS PROCEDURE IS VITAL TO 


THE OPERATION OF ''TIE PROGRAM 
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'PACK' IS WRITTEN IN UNBLOCKED BCD CARD IMAGES, EVEN PARITY, 800 BPl 
DENSITY, ON A T/\PE LABELED 'PACK TEMP. PROGRAM'. IT SHOULD BE COPIED TO 
DRUM AND ANY REVISIONS HANDLED ON THE DRUM COPY, LEAVING THE TAPE COPY IN 
THE ORIGINAL FORM. IN THE PRESENT FORM, THE PROGRAM REQUIRES OVER 6500 
LOCATIONS FOR INSTRUCTIONS AND OVER 6600 FOR DATA (INCLUDING ALL LIBRARY 
ROUTINES) . 

OUTPUT TAPE FORMAT 

pack' WRITES AN OUTPUT TAPE WITH 755 WORDS (36 BITS) PER RECORD, AT 
800 BPI DENSIT-Y AND ODD PARII’Y. EACH RECORD THEREFORE OCCUPIES 5.6625 
INCHES OF TAPE, PLUS 0.75 INCHES FOR IRG, FOR A TOTAL OF 6.4125 INCHES 
PER RECORD. EACH INPUT TAPE HAS ABOUT 2400 RECORDS WITH 8.98246 SCANS 
PER RECORD, SO FIVE INPUT TAPES HAVE ABOUT 107,790 SCA-NS. EACH OUTPUT 
RECORD HAS 25 SCDVNS , SO THE OUTPUT TAl’E WILL CONTAIN ABOUT 4315 RECORDS, 
USING APPROXIMATELY 2,306 FEET OF TAPE. 

THE FIRST 5 WORDS OF EACH RECORD CONTAIN CONSTANTS AND FACTORS: 


RECORD 

WORD 

1 

BITS 

1 - 

9 

REEL NUMBER 




BITS 

10 - 

18 

TA’E TYPE 




BITS 

19 - 

36 

RECORD NUMBER 

RECORD 

WORD 

2 

BITS 

1 - 

9 

YEAR ( 9 OR 0) 




BITS 

10 - 

18 

DAY 




BITS 

19 - 

36 

OCTAL 300000 

RECORD 

WORD 

3 

BITS 

1 - 

18 

FACTOR 3187 




BITS 

19 - 

36 

FACTOR 52 


RECORD WORD 4: SAME AS SCAN WORD 1 FOR THE FIRST SCAN IN THE RECORD. 
RECORD WORD 5: SAME AS SCAN WORD 1 FOR THE LAST SCAN IN THE RECORD. 







FOLLOWING THE FIVE HEADER WORDS ARE 25 SC/VNS OF DATA, ALL CODED THE 


SAME, 

WITH 

30 

WORDS 1 

PER 

SCAN. 


SCAN 

WORD 

1 

BITS 

1 

- 

6 

FI.AG 




BITS 

7 

- 

18 

DAY 




BITS 

19 

- 

24 

HOUR 




BITS 

25 

- 

30 

MINUTES 




BITS 

31 

- 

36 

SECONDS 

SCAN 

WORD 

2 

BITS 

1 

- 

18 

FIRST TEMPERATURE 




BITS 

19 

- 

36 

SECOND TEMPERATURE 

SCAN 

WORD 

22 

BITS 

1 

- 

18 

41 ST. TEMPER.\TURE 




BITS 

19 

- 

36 

42 ND. TEMPER.\TURE 

SC/\N 

WORD 

23 

BITS 

1 

- 

18 

FIRST WIND 




BITS 

19 

- 

36 

SECOND WIND 

SC,VN 

WORD 

30 

BITS 

1 

- 

18 

15 TH WIND 




BITS 

19 

_ 

36 

16 TH WIND 


(NOTE: BITS ARE NUMBERED FROM LEFT TO RIGHT) 

THE FLAG IS SET TO "1" FOR A NORMAL SCAN; IF A SCAN IS MISSING OR IN ERROR, 
THE FLAG IS SET TO "2", THE PROPER TIME FILLED IN ,\ND ALL DATA POINTS SET 
TO ZERO. IF THE IJiST INPUT SCAN FROM THE FIFTH INPUT TAPE DOES N'OT RESULT 
IN A COMPLETE OUTPUT RECORD, THE PROGRAM ZERO FILLS THE Rj;ST OF THE RECORD, 


INCLUDING FLAG ,VND TIMli BITS 


r 


If) 


WHEN THE FIVE INPUT T,\PES ,\RE FINISHED, THE PR(X;R.VM WRITES AN EOF, 
THEN WRITES A TR/MLER; TRAILER WORD 1: SAME AS RECORD WORD 2 FOR HiE 
FIRST RECORD ON THE TAPE. 

TR^ULER WORD 2: SAME AS SCAN WORD ] FOR THE FIRST SCAN ON THE TAPE. 

TRAILER WORD 3: SAME AS RECORD WORD 2 FOR THE LAST RECORD ON THE TAPE. 

TRAILER WORD 4; SAME AS SCAN WORD 1 FOR THE LAST SCAN ON THE TAPE. 

TRAILER WORD 5: SAME AS RECORD WORD 1 EXCEPT WITH TOTAL RECORD COUNT. 

THE PRO(;rAM then writes TWO MORE EOF'S. 

SAMPE RECORD: 


WORD 1 


5 

7 

27 

35 

751 

753 

755 


9 

REEL 


9 

TYPE 


18 

RECORD NO. 


9 : 9 ; 18 j 

jYEAR ; DAY joCTAL 300,000 j 


I 18 
DECIMiU. 3187 


18 

DECIMAL 52 


SAME AS WORD 6 



6 

8 

2B 

36 

752 

754 




EACH KEtlORl) STARTS WITH THE FIVE WORD HE;\DER; THE TRAILER RECORD IS WK I 
ONLY ONCE PER TAPE, FOLLOWINC. AN END-OF-FILE MARK. 

TR.\ILER RECORD 


EOF 


• 9 t 9 1 18 j 

1 1 VR ' DAY 'OCTAL 300,000 | 

6 I 12 ' 6 • 6 ' 6 ; 

Fl.\c;' DAY 1 HR 1 MIN 1 SEC j 

[ FROM FIRS! RECORD ON THIS TNU’E ] 

9 I 9 1 18 I 

} I YR ! DAY loCTAL 300,000 

b ! 12 I 6 I 6 ; b 

FLACI DAY 1 HR MIN ! SEC j 

1 FROM LAST RECORD ON THIS TAPE ] 

9 9 18 

5 1 REEL TA’PE TOT,\L RECORD COUNT 

ZERO FILL OR TR^YSH ^ 


7 ! ZERO fill or trash 

1 . . . _J 

ZERO FILL OR TR,\SH j 

75 3 j ZERO fill' dlTmYsT! ' 

ZERO FILL OR TR^YSll 
_1 

755 ZERO fill OR ~f RASH 

1 


“'eof' 

EOF 

TO CONVERT THE DATA TO USABLE FORM, SUBTRACT llCTAL 300,000 FK"" 


rEMl’ER.\TURE POINT ,VND MULTIPl.Y BY O.OI. EACH WIND POINT MUST BE 
A FU'ATINC POINT NUMBER, DIVIDED BY 31.87, AND THE RESULT ADDED TO 0.3 
OBTAIN MILES PER HOUR. AS M.\NY OPERATIONS AS POSSIBLE SHOULD BE D.'-;! 
niE fixed POINT REPRESENTATION SINCE DIVIDE IS ,VN EXPENSIVE OPl.RAIl.'N. 


L IL 


1 « 


NOTE ON READING TAPES ON GDC 6600 
EACH RECORD CONTAINS 453 WORDS (60 BIT) , WITH 18 WORDS PER 
RECORD SHOULD BE BUFFERED INTO ,VN UNPUT ARRAY, THEN TRAI.'SFER TO 
PROGRAM. SINCE FIVE '1108' WORDS EQUAI. THREE '6600' WORDS, THE 
CAN WORK IN UNITS OF THREE WORDS. 


SC/VN. THE 
AN ASSEMBLY 
PROGR.\M 



19 


In trying to iraplenent this program on the National Bureau of Standards' 
UNIVAC 1108 under LX£C 2, several problems were encountered. The program 
would have tied up six tape drive units during execution but only two at a 
time are actually used, so the systems staff at NBS suggested using one 
multi-reel file for input ratlier tlian rive. After this cliange was made, 
the program would not execute properly. The problem was finally traced to 
the 'N'TRiVi' tape swapping routine which works properly under KXhC 8 but not 
in EXEC 2. The systems staff at NBS was very helpful but none of their 
suggestions resolved the difficulty. Therefore, tlie only alternative was 
to vjrite the tape swapping section directly in assembly language. However, 
this would have tied the program to EXEC 2 since the tape routines di^'fer 
between the two versions of EXEC. Since N’BS was in the process of converting 
to EXi-C 8, a more logical procedure was to use EXEC 8 and not try to make 
the progr;im compatible with EXEC 2. Due to various logistics problems 
involved with D-tape data reduction and processing off campus, a decision 
was made during the fourtli quarter to perform this task at the University 
of .'laryland Computer Center rather than at the National Bureau of Standards. 

Another problem came from the 14U1 tape copies of the original data 
tapes. For some reason, the copied tapes liad dropped bits and some evidence 
of "print througli" or bleeding of bits was indicated. Tlie 'NTR/\N' routine 
was not able to handle these problems economically, so a liardware-level routine 
was inserted for the'HTRAN' read routine and direct assembly routines used 
for tape swapping and output. 

A section was added to the input routine to make possible direct 
reading of fairly good tapes witliout going tlirougli tiu' 1401 copv process. 
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Tills routine will work wltli most tapes as long as there are no end-of- 
files except at the beginning and end of the tape. However, the initial 
tests with poor tapes showed that the 1401 could recover more of the data 
than the 1108. 

A program was written to process the output tapes from 'PACK'. The 
program computes half-hour average temperatures and wind speeds centered at 
the half-hour and hour. Since one tape provides about two weeks of data, 
longer averages are also possible without tape swapping routines. The 
analysis in Section IV makes use of such 30-minute averaged profiles except 
for certain of the spectral analyses. 

Considerable progress was made during this project in converting the 
raw TREMD data tapes into usable Information. This progress was not without 
its problems and delays, as noted before. The efficiency of data reduction 
increased witli the transfer of all processing to the University Computer 
Center althougli the unanticipated extra expense of preliminary IBM 1401 
copying and Univac 1108 repacking was borne by this contract account. This 
ultimately made it necessary to terminate these preprocessing activities 
before all the raw D-tapes were repacked. In retrospect, however, we are 
quite pleased to have salvaged as mucli data as we have. 

Table 1 provides a documentation of the D-tape data processing. Tlie 
set of remaining raw input tapes is "rough" enough to cause minor problems 
in any future attempts to complete the preprocessing of the TIUCMD project 
data record. For example, certain raw I)-tapes (e.g.: //322,388) are missing 
so that there will be gaps in the data coveragi' unless they are found. Other 
input tapes (e.g.: #340) yielded far less data than expected after iirocessing. 




TABLE 1. SUMMARY OF PROCESSED D-TAPE 

DATA (Both 7-track 

and 9-track output 

capes) 


<eel 

.umoer 

Input Tape: 

s 

Dates of 

Data Coverage 

Input Records 

Output Records 

Total Scans 

Bad Scans 

y Bad 
Scans 

J40 

237, 238, 
241, 243 

240, 

(9/25/69) - 
(10/8'69) 

12,112 

4314 

107,850 

1073 

1.0 

J 73 

273, 274 , 
277, 278 

275, 

(12/1/69) - 
(12/12/69) 

10,871 

3879 

96,975 

110 

O.Il 

1^7 

283, 286, 
288, 280 

287 . 

(12/26/69) - 
(1 '8. 70) 

10,717 

3850 

96,250 

8 

Q.Ol 

39J 

291, 292, 
293, 29n 

293, 

(1'S'70) - 
(1/19.'70) 

10,617 

3814 

95,350 

10 

0.01 

Wl 

298, 299, 
302, 304 

301 , 

(1/19/70) - 
(2/13/70) 

10,938 

3946 

98,650 

409 

0.41 

308 

305, 306, 
309, 311 

iOH 

(2 '13/70) - 
(2/28, 70) 

12,006 

4306 

107,650 

79 

■1.07 

ilh 

312, 3U, 
317, 318 

iW), 

(2/28/70) - 
(3M4/703 

11,962 

4240 

106,000 

273 

0. Jb 

331 

328, 330, 
333, 334 

)3l. 

(4/6/70) - 
(,./20/70) 

12,048 

431 1 

107,775 

134 

0. IJ 

33S 

333, 337, 
34 0, 341 

333, 

(4/20/70) - 
(5/3/70) 

1 1 .092 

3862 

96,550 

617 

O.'M 

334 

351, 352, 
355, 33i, 

334 

(5/24/7U) - 
(6/6/70) 

12,127 

3 3 55 

83, .875 

16,121 

]‘i. J 

3(.3 

363, 363, 
367, 169 

luu. 

(6/19/70) - 
(7/2/70) 

11,750 

4096 

102,400 

1,271 

1.24 

374 

371, 372, 
175, 376 

374 . 

(7, 2 '70) - 
(7/24/70) 

10,749 

3862 

96,550 

38 

0.0.', 

3.30 

377. 378, 
381, 383 

380, 

(7/24/70) - 
(8/11/70) 

12,000 

4287 

107,175 

7 

•7 

3‘35 

393, 39-., 
397, 398 

395 , 

(8/29/70) - 
(9/13/70) 

11 ,981 

3972 

99,300 

21,729 

21,8s 

i... 

/' J 



1o:,97j 

3b . 094 

1 .402,350 

41,872 

2 . OT 


il LURKS; 

Input 

I n pu t 
Input 

tape 278 Is short, 
tape 295 Is short, 
tape 298 covers the 

TOT,\LS 

period (1/19/70) - (1/21/70) then 





.skips to (2/5/70) 

Input tapes 322 and 388 are mlsalnK* 

Input tape 323 has nothing recorded on it. 

Input tape 340 is short. 

Input tapes 358, 362', 362" have special problems. 

Input tapes 384, 385. 389 liave "noise" problems related to 
simultaneous recording of A-tape data. 

Input tape 393 yielded only 3 good scans. 

Input tape 394 yielded about 13,000 of 19,000 scans. 


Sucii tapes had to be recopied and reprocessed a second time to obtain more 
nearly the full ;imount of data. Further, at least one raw input tape 
Otili) was found with no data on it at all. Some input tapes (e.g.: '*384, 
385, 389) presented repacking problems due to "noise" v;hich was generated 
by simultaneous A-tape recording on site. On the more encouraging side, 
some raw input tapes were repacked without the preliminary IB>! 1401 copy 
w'ork. Of the 40 input tapes that have been fully processed thus far, 

810,075 scans have been recovered of which only about 2-1/4"' are bad or un- 
usable scans. Overall, our experience v;itii tlie data reduction indicates 
that a large faction of the TREND project micro-meteorological tower data 
can be salvaged from t'ne raw data tapes by the processing routines developed 
by our research group provided proper care and time is given to the effort. 

It should be mentioned that final output tapes are \.'ritten in botli 7- 
track and 9-track versions. Since it is difficult to make 7-track copies 
of 9-track tapes, tliis dual output was done in order to preserve the unique 
TREND data as a universal data source for future analysis and research. 
Hence, any future requests for copies of the TREND micro-meteorological data 
can be handled using output tapes of the s^tme format as requested. 

It was virtuallv impossible to prepare analyzable data tapes from 
the analogue (A) tapes partly due to extra financial squeeze and partly 
due to insufficient information regarding tlie original preparation of 



A- tapes . 
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Certain radiation, dew point temperature and wind direction information 
for limited periods have been keypunched from data printouts made from 
some of the A-tapes during the original TREND field experiments. These 
data are summarized later in this report. In addition, some data from the 
various ground stations in the TREND experimental area were made available 
to us by USAETL and these serve as the principle source of information 
for the climatological summary in the next section. 


SECTION III - CLIMATOLOGICAL CONDITIONS IN THE SAKAERAT FOREST 

1. CLI.'LU'OLOGICAL SUMMARY OF THE EXPERIMENTAL AREA 

To place the micrometeorological studies of the forest into proper clima- 
tological context, data from the principal surface station were summarized 
over tl\e nearly three years of record to construct a climatology of the area. 
Computer printouts of (lourly values of temperature, relative humidity, pre- 
cipitation, pan evaporation and hours of sunshine were available from tl\e 
original field experiment. The data were averaged over each month of the 
record from November, 1967, to September, 1970. Using the average air tem- 
perature and relative humidity values for each month, values of "effective 
temperature" were obtained in conventional fashion (Landsberg, 1970). Basi- 
cally, the effective temperature takes into account tlie comtiined effect of 
temperature and humidity on a person's comfort. Effective temperature is 


! 

Jit 

I 

defined as a sensation in which, for a given air temperature and humidity, 
the state of comfort is equal to that experienced for an environment at a 
lower temperature with saturated conditions. That is, air would have an 
effective temperature of 20°C if it feels the s;ime as air at 20°C and 100% 
relative liumidity. For example, this would be the case for air at 22.5°C 
and 50% relative humidity or at 25°C and 10% relative humidity. Averages of 
hourly values of effective temperature were compared with values computed 
from daily averages of temperature and relative humidity for three test days — 

December 14, 1967; April 19, 1970; August 13, 1970 — to see if one could ob- 
tain reasonably close values with the lesser computation frequency.* The mean 
daily effective temperature for tiiese three days were, respectively, 14. 3°C, 

27.1°C and 24.7°C. Corresponding values of effective temperature computed 
from the daily mean temperature and relative humidity were, respectively, 

14.6°C, 27.4°C and 25.1°C. These are tolerably representative values so the 
sumraaries incorporate values computed on this latter basis. 

Precipitable water, U, in the atmospheric column over Sakaerat Experi- 
ment Station were computed using the metl of Smith (1966) with empirical 
(■) values taken from Lettau (1970). These were computed using monthly 
average values of temperature and relative, humidity. Values of tlie storage 
of water vapor in the atmosplieric column were computed by taking the diffe- 

I 

rences of precipitable water from month to month, AU. The advection. A, I 




* See Section II 1.6 for furtlier ilocument at ion . 
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of water vapor into the atmospheric column was computed as a residual in 
the balance equation. 

Storage = Evaporation - Precipitation + Advection 

AU = E - P + A 

Table 2 gives values of climatological elements for each month of the 
record at Sakaerat Forest. Figure 4 shows a representative annual variation 
of each element determined by further averaging overall corresponding 
months in the record. The figures on the number of hours of daylight were 
taken from the Smithsonian Meteorological Tables on Duration of Daylight 
for the 15th of each month. 

The cool, dry northeast monsoon, associated with high pressure over 
the Asian continent, is evident in each year of record in Table 2 during 
November to February as well as in the annual summary. During these periods, 
the temperatures and relative humidities are comparatively low. The tempera- 
ture values indicate that the northeast monsoon was less vigorous during 
1968-69 than the other two years in the record. The average temperature 
during these four months of 1968-69 was 25.4°C as compared with 21.7°C for 
1967-68 and 22.2°C for 1969-70. 

The representative annual cycles of temperature and effective tempera- 
ture show maximum values (27.5°C, 24.9°C) during the transitional month of 
March and minimiun values (21.5°C, 20.0°C) in December. The annual average 
temperature at Sakaerat was 25.0°C and annual average effective temperature 
was 23.5°C, with only a 5-6C° variation in either during the year. 

* These data are plotted in Figure 9. 




TABLE 2. Monthly Values of Climatological Elements at 
Sakaerat Experiment Station, Thailand. 


Month 

T 

CC) 

R.H. 

(%) 

T 

eff 

(°C) 

Hours 

Sun- 

shine 

Hours 

Day- 

lip.ht 

Precipitation 

(mm) 

Total Daily 

Evaporation 

(mm) 

I’o t a 1 Dally 

11 

(mm) 

AU/At 

(mm) 

A 

(mm) 

11-67 

22.74 

85.0 

21.8 

6.4 

11.43 

92.4 

3.1 

- 

- 

37.7 


- 

- 

12-67 

19.13 

77.0 

18.2 

7.6 

11.25 

00.0 

0.0 

- 

- 

27.0 

- 

10.7 

- 

1-68 

21.12 

68.8 

19.9 

7.7 

11.33 

00.0 

0.0 

- 

- 

27.4 

+ 

0.4 

- 

2-68 

22.82 

67.4 

22.4 

7.8 

11.67 

2.1 

<0.1 

- 

- 

33.6 

+ 

6.2 

- 

3-68 

26.36 

70.8 

24.3 

7.7 

12.03 

23.2 

0.75 

- 

- 

38.4 


4.8 

- 

A-68 

25.38 

78.6 

23.9 

6.3 

12.47 

138.7 

4.6 

- 

- 

40.2 

+ 

1.8 

- 

3-68 

26.09 

82.8 

24.8 

6.3 

12.83 

227.2 

7.3 

- 

- 

44.2 


4.0 

- 

6-68 

26.28 

82.4 

25.0 

4.4 

13.00 

81.0 

2.5 

- 

- 

48.4 


4.4 

- 

7-68 

26.58 

78.8 

25.3 

5.5 

12.93 

32.5 

0.9 

147.56 

4.76 

47.1 

- 

1.3 

-116.4 

8-68 

26.59 

77.3 

24.9 

4.9 

12.63 

50.9 

1.5 

137.70 

4.59 

46.1 

- 

1.0 

- 87.8 

9-68 

25.43 

78.5 

24.4 

5.2 

12.22 

138.3 

4.6 

111.00 

3.70 

41.1 

- 

5.0 

22.3 

10-68 

23.47 

85.0 

22.5 

6.3 

11.82 

89.4 

2.9 

95.48 

3.08 

39.7 

- 

1.4 

- 7.5 

11-68 

24.37 

79.2 

23.0 

7.9 

11.43 

4.9 

0.1 

125.70 

4.19 

39.1 

- 

0.6 

-121.4 

12-68 

24.95 

67.3 

22.8 

8.1 

11.25 

0.0 

0.0 

145.08 

4.68 

33.9 

- 

5.2 

-150.3 

1-69 

25.12 

68.9 

23.0 

6.5 

11.33 

15.3 

0.5 

122.45 

3.95 

34.9 


1.0 

-106.2 

2-69 

27.24 

60.3 

23.4 

8.1 

11.67 

0.0 

0.0 

183.40 

6.55 

34.6 

- 

0.3 

-183. 7 

3-69 

27.58 

65.5 

24.9 

6.9 

12.03 

152.9 

5.0 

175.15 

5.65 

38.2 


3.6 

- 18.7 

4-69 

27.53 

68.3 

25.0 

7.5 

12.47 

74.1 

2.5 

187.20 

6.24 

39.6 


1.4 

-111.7 

5-69 

27.02 

76.2 

25.1 

5.6 

12.83 

110.1 

3.7 

150.04 

4.84 

42.9 


3.3 

- 36.6 

6-69 

26.53 

82.1 

25.2 

4.6 

13.00 

146.0 

4.9 

137.70 

4.59 

48.7 


5.8 

14,1 

7-69 

25.71 

82.6 

24.5 

4.1 

12.93 

110.3 

3.5 

122.14 

3.94 

46.8 

- 

1.9 

- 13.7 

8-69 

25.62 

80.9 

24.3 

5.4 

12.63 

169.6 

5.5 

107.88 

3.48 

45.5 

- 

1.3 

60.4 

9-69 

25.02 

85.0 

24.0 

4.2 

12.22 

370.1 

12.3 

72.60 

2.42 

43.8 

- 

1.7 

295,8 

10-69 

24.48 

85.0 

23.5 

5.3 

11.82 

151.4 

4.9 

84.94 

2.74 

42.2 

- 

1.6 

64,9 

11-69 

21.59 

81.6 

20.6 

6.2 

11.43 

57.4 

1.9 

76.80 

2.56 

34.0 

- 

8.2 

- 27.6 

12-69 

20.36 

70.3 

19.0 

7.4 

11.25 

missing 

121.83 

3.93 

26.6 

- 

7.4 

- 

1-70 

21.70 

69.5 

20.1 

7.4 

11.33 

18.8 

0.6 

127.41 

4.11 

28.7 


2.1 

-106.5 

2-70 

25.29 

63.3 

22.8 

8.2 

11.67 

10.3 

0.4 

172.20 

6.15 

32.5 


3.8 

-158,1 

3-70 

28.45 

65.5 

25.6 

7.2 

12.03 

71.5 

2.3 

213.28 

6.88 

40.2 


7.7 

-134.1 

4-70 

26.89 

73.3 

24.8 

7.1 

12.47 

136.3 

4.6 

149.40 

4.98 

41.0 


0.8 

- 12.3 

5-70 

27.21 

76.9 

25.4 

6.1 

12.83 

177.4 

5.7 

141.40 

4.55 

43.8 


2.8 

38,8 

6-70 

26.62 

79.5 

25.1 

5.1 

13.00 

139.8 

4.7 

113.40 

3.78 

47.4 


3.6 

30.0 

7-70 

25.98 

79.6 

24.5 

4.5 

12.93 

90.7 

2.9 

129.27 

4.17 

45.9 

- 

1.5 

- 40.1 

8-70 

25.26 

81.5 

24.0 

4.9 

12.63 

146.0 

4.7 

100.75 

3.25 

45.1 

- 

0.8 

44.5 

9-70 

25.11 

83.1 

24.0 

5.0 

12.22 

248 . 3 

9.2 

95.10 

3.17 

42.8 


2.3 

155.5 
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The relative humidity generally falls during the monsoonal period 
reaching a minimum value in February. The period minimum value was slightly 
more than 60% during February, 1969, the warmest February in the record. 

The relative humidity is generally largest in October, with an annual average 
value of 85%, and smallest in February with annual average value of about 
64%. The largest difference between air temperature and effective tempera- 
ture occurred during February, 1969 when the relative humidity reached its 
period low value. The average air temperature during that month was 27.2°C 
while the effective temperature in the relatively cool dry monsoonal air 
was only 23.4°C, almost 4C° cooler. In fact, the effective temperatures 
during the 1968-69 northeast monsoon and continuing through the first six 
months of 1969 were considerably cooler than air temperatures due to a long 
spell of relatively dry air associated with that season's monsoon. 

The northeast dry monsoon occurs during a period of relative minimum 
day length yet the sunshine recorders measured relative maximum amounts of 
sunshine during this period which indicates smaller amounts of cloud cover. 

The precipitation traces are also small averaging generally less than 1 mm/day 
during the northeast monsoon. 

The warm and moist southwest monsoon generally occurs during the period 
April to September corresponding to the months of maximum dayliglit. The 
temperature, precipitation and relative humidity are generally high during 
this period. The total rainfall during these six months was 980 mm for 1969, 
939 mm for 1970 and a relatively smaller amount of 669 mm during 19b8. Con- 
tributing to the large total for 1969 was an average daily rainfall during 
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September of that year of nearly 12 mm/day. The double maxima in precipita- 
tion during the southwest monsoon, shown in Figure A, is also evident in a 
siunmary of Bangkok climatological data for the period 1943-1952 from the 
noyal Thai Navy which was prepared by Lettau (1970). This will be discussed 
further in the next section. 

The relative humidities are high during these monsoonal periods yet the 
effective temperatures are lower than the air temperatures by a greater amount 
than during the dry northeast monsoon where the air temperatures are much lower. 
The humidity generally increases during the southeast monsoon period achieving 
its tnaxinuim toward tlie end of the period, althougli there was an intrusion of 
drier air during July-September of 1968 which modified this trend. 

The higher degree of cloudiness associated with the soutliwest monsoon 
is reflected in the sunshine recordings where relatively smaller amounts of 
sunshine occur each day than during the northeast monsoonal period. In 
fact, the monsoon cycle is well illustrated by the sunshine recording whiclt 
Is roughly 180° out of phase with the cycle of daylight duration. 

The precipltable water record tends to parallel the relative humidity 
record over the whole period. The effect of temperature on the precipitablc 
water is best sliown for September-Novenber , 1968 where the lowering air 
temperature maintains a lowering precipltable water value in the face of a 
large peak in relative humidity values. The evaporation tends to respond 
to the solar radiation available at tlie surface but is also controlled by 
the ambient relative humidity and air temperature. The evaporation is 
generally largest wlien air temperatures are higlt and ^imbient hiunidity is 
low, and generally smallest iluring periods of Iiigher ambient luunidity, 
falling temperatures and fewer liours of sunsliine. Similar tenper.ature effect-. 
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occur during the same months of 1969 and 1970 and seem somewhat characteristic 
of the monsoonal climate. 

The water vapor storage term is generally small throughout the period 
with values within + 10 mm/month and is not shown in Figure 4. The advec- 
tion term shows the large monsoonal influx of moist air feeding the rainfall 
systems of May and September with outfluxes of moist air from the region 
during periods of evaporation excess over precipitation during the northeast 
monsoon. Average rainfall during May and September is 172 inni/month and 
252 mm/month respectively separated by three months of lesser rainfall 
with local minimum during July of 78 mm/montli. The precipitation drops to 
near zero during the northeast monsoon season of December, January and 
February. For the rainfall peak of May, about 90% (14b mm) is contri- 
buted by evaporation, 17% (29 mm) is contributed bv advection with 3. 3 mm 

added to storage. During September, 62% (156 mm) is contributed by 

advection, 37% (93 mm) is contributed by evaporation and 1.2% (3.0 mm) 

is contributed by storage. During the relative minimum of July, there is a 
fairly vigorous advection of about 43% (57 mm) of the evaporated water 

out of the region leaving only about 78 mm to rain out. During the dry 
montlis of December, January and February nearly all tlie evaporated water 
is advected out of the column with little precipitation and only minor 
clianges in storage. 

2. DIURNAL VARIATION OF CONDITIOIIS AT S/iKAFRAT FOREST 

Data from four different months of the year were analyzed for their 
diurnal variations. Figures 5a, b, c, siiow the corresponding values of 
temperature, relative hiunidlty and precipitation on an hourly basis for the 




The Day 


months of March, June, September and December. Directing attention first 



to the temperature. Figure 5a shows its diurnal variation for each of tliesc 
months. Temperature maxima generally occur around 1300 LST-IAOO LST. i'iu 
diurnal variation for the month of September shows a maximum somewhat ear- 
lier than this. Referring for the moment to the precipitation record in 
Figure 5c, it is seen that the rainfall during September has a maximum at 
1400 LST with much lower values during the preceding two hours but maintai:'. 
high values throughout the afternoon. This characteristic of the diurnal 
rainfall cycle effectively shifts the temperature maximum to an earlier i.; ■. 
by cooling the air during those rainy hours when the temperature would 
ordinarily be expected to achieve its maximum value if skys were clear. 

As was apparent in tlie monthly temperature summary, the temperature i 
highest during the transi*'ional montli of March wiien there are more hours . ■ 
sunshine near the surface. Temperatures are lowest during December when 
the solar cycle yields the shortest daylight period. For eacli of tiie mom 
shown, temperatures liave their minimum values around 0600. 

The relative humidity during tlie southwest monsoon months is considi i 
ably higher than during the transitional month of March and tlie month ot 
December during the dry northeast monsoon. During die months of March aiii: 
December, relative humidity reaches its minimimi value around 1400-1500 l.Si, 
the minima during June and September occurring earlier, alunit noon. In ■ .■ 
montli, the minima in relative humidity tend to occur in association with 
nuiximum temperatures when saturated mixing ratio values would achieve their 
maxima also. Relative humidity is generally largi'st in the predawn lieut 
coinciding with the time of minimum temperature and saturated mixing i ii . 
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The precipitation during December is virtually zero at all hours. 
Afternoon maxima are evident in each of the other months shown with September 
and Jiuie showing tlie largest values. There appears to be something of a short 
periodic cycle in the afternoon rainfall, although we did not statistically 
test for tills. During the rainy season, September shows three local peaks 
in the afternoon occurring at 1200, 1400, and 1600 and June shows local 
peaks at 1300, IbOO and 1900. This suggests a 2-3 hour cycle in the after- 
noon monsoonal convection. Each of the three months of significant rain also 
show nocturnal maxima; a peak occurring in the early evening between 2000 and 
2400, depending on the month, and an early morning maximum a few hours before 
dawTi. All of the months show relatively small amounts of rainfall between 
dawn and noon. 

3. INTERDIURMAl. TEMPERv\TURE VARIATIONS 

The change of temperature from day to day is a good measure of climatic 
conditions. Ft gives a good index of restlessness of tlie local atmospFiere. 

A frequency distribution of this element for Sakaerat Forest is sliown in 
Table 3. The table sliows that the frequency of rises and falls of temperature 
from day to day is nearly symmetrically distribvited around zero. The lack 
of prominent synoptic disturbances is exetiqilified by tlie fact that 90/i of 
the interdiurnal clianges are less than 2°('.. Tliese are essentially imper- 
ceptible and characterize a monotonous bioclimate witliout stimuli. In the 


whole period of record tliere were only 8 cases wlien rises and falls between 
days equalled or exceeded U°C. Six of these were falls and are minor cold 
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outbreaks which once in a while even befall these low latitudes. By middle 
latitude standards one would hardly notice them, considering the fact that 
the mean interdiurnal temperature variations on the U. S. Great Plains In 
November to April is about 5°C and that values of 3°C are very common 
(Landsberg, 1966). 

The monotony of the climate is also borne out by an iteration analysis 
(see Yule, 1950). This consists of a check on the length of period for which 
the same sign of interdiurnal variation persists. Table 4 gives the number 
of cases for rises and falls of temperature of varying length in days, due 
allowance being made for missing data. There is a tendency for rise periods 
to be somewhat longer than the falls, indicating that any cooling trend is 
more apt to be short lived than warming trends. 

The iteration test yields the results of Table 5. Even without any 
further test it is immediately noted that the quick iterations from falls 
to rises are considerably less than statistically expected and the long 
intervals without iteration in sign of temperature trend quite inflated. 

In fact, a chi-square test indicates a probability of 0.005 that the ob- 
served distribution could be produced by chance and deviate this much from 
the expected. The statistical test confirms formally the earlier statement 
that day-to-day changes tend to persist in sign and also that this is more 
pronounced for rises than falls of temperature. If one compares this type 
distribution with conditions in higher latitudes, especially in winter, one 
encounters there a great many cases in the 1-day category with sign changes 
in the short duration category overwhelmingly common. Instead of the 
tendency for persistence of trend, as at the tropical station, there is a 
distinct pattern of oscillation. Tlie period of observation at Sakaerat 
Forest is too short to attempt a seasonal breakdown. Similarly the data am 
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Tabic 3. Percent Fretiuency of Interdiurnal Temperature Change at 

Sakaerat Forest. 


“C 0 0-0.9 1.0-1. 9 2.0-2. 9 3. 0-3. 9 

L rises 'i 32 13 4 1 1 

.i lali^ 0 2C 13 3 10 


(976 cases) 
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33 
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0 
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0 
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33 

23 
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1 
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1 
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fable 5. Iteration Test 
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not adequate to perform an harmonic analysis in the search for possible 
rhythms in weather patterns. 

4. COMPARISON BETWEEN SAKAERAT FOREST AND OTHER SOUTHEAST ASIAN LOCATIONS 

Data presented by Anstey (1966) for several locations in Southeast Asia 
were compared with that at Sakaerat Forest to determine certain geographic 
differences as well as differences due to forest exposure of instrumentation. 
Table 6 shows the geographic configuration of all stations and the annual 
total rainfall of each. Table 7 and Figure 6 show the comparison of basic 
elements at four of these stations which are at about the same latitude. 

The temperature record at Sakaerat Forest closely parallels that for 
Bangkok, which is about 190 km away from the TREND site. There is a fairly 
consistent difference in temperature between these two locations which 
averages 2.5°C through the year, Bangkok showing the warmer temperatures. 

This difference is due largely to differences in elevation, instrument ex- 
posure between the two locations and the city influence. Subtracting this 
annual average difference of 2.5°C from Bangkok values leave remainders 
which are within 1.1°C of the Sakaerat readings for virtually every month 
of the year. This is so even though the data for Bangkok are from different 
years than those for Sakaerat Forest. The main differences in the annual 
variation in temperature from the Bangkok summary and the Sakaerat Forest 
data occurs during April and May, at the beginning of the moist, southeast 
monsoon season. The correlation of temperature data at Sakaerat Forest witli 
those at Nhatrang and Tavoy are generally worse, especially with Tavoy sliowin 
a much different annual cvcle of temperatures. 
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TABLE b. GEO(;iL\PHICAL CONFIGURATION OF SOUTHEAST ASIAN LOCATIONS 
RELEVANT TO THE CLIMATOLOGICAL SUI’ ^IA' IJIS ''F TABLE 7, FIGURES b, 7, 8. 


Annual loiai 


Location 

Elevation (M) 

Latitude (N) 

Longitude (E) 

Rainfall 

Bangkok (Thai land 

2 

13°44' 

100°30' 

1400 

Ciiapa (N. VietNaiu) 

lb41 

22° 22 ' 

103°52' 

2860 

Dalat (N. VletlLun) 

1501 

Il°57' 

108°2b' 

1803 

Hanoi (N. VietNam) 

7 

2i°03' 

105°52’ 

1808 

Kota Bharu (Malaya) 

b 

b°08' 

102°15' 

3147 

Luang Prabang (Laos) 

340 

1V°53' 

102°08' 

1308 

Mandalay (Burma) 

7b 

21°59' 

9b°06’ 

84b 

Nliatrang (S. VietNam) 

b 

12°15' 

109° 12' 

1440 

Rangoon (Burma) 

b 

lb°47' 

9b"13' 

250 7 

S.ilaoa Is. I'iftN.im) 

11 

10°47’ 

10b°40' 

19 74 

Sakaerat Forest (Thai land) 

550 

14°31' 

101°55’ 

113 3 

Tavoy (Burma) 

h 

I4°05’ 

98°12' 

5291 
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The record of relative humidity at Sakae.rat Forest most closely re- 
sembles that of Bangkok but in this case, the discrepancies are much larger. 
Over the year, the relative humidity at Sakaerat Forest averages 9 percentage 
units higher than Bangkok but the month to month differences vary from 1 
percentage unit in January to 14 percentage units in May, August and November. 
The differences in climatic regimes at the other two locations is again 
in evidence in the relative humidity records with Nhatrang exhibiting a 
somewhat similar regime during the southwest monsoon season but departing 
considerably during the northeast monsoon and spring transition. Tavoy, 
once again, shows a completely different type of annual variation. 

The precipitation records of Bangkok and Sakaerat Forest show similar 
regimes during the entire year with monsoonal peaks during May and September, 
with an interceding relative minimum, and a comparatively dry winter. 

Averaged over the year, monthly precipitation at Bangkok is about 22 ram 
higher than at Sakaerat Forest, however, the range of monthly differences 
for individual months throughout the year is from 85 mm to -58 mm. The 
precipitation at Nhatrang shows a large monsoonal surge from August to 
December peaking in November and a substantially smaller, early season 
peak in phase with that in Thailand. The precipitation at Tavoy, Burma, 
which is not shown in Figure 6, exhibits only a single surge during the 
period April to November with a maximum monthly rainfall of 1242 mm occurring 
in July. The comparatively large values at Tavoy are presumably dtie to the 
lack of mountains or coi.tlnental area on its winci^ard side during the 


southwest monsoon. 



Since some but not all stations in Anstey's summary exhibit double 
maxima in precipitation, a further comparison of the rainfall regimes of 
, these locations was made in search of a consistent pattern. Figures 7a-e 

f show a normalized [percent of annual total] annual variation of precipitation 

for selected stations chosen to illustrate the various monsoonal regimes. 

The data have also been normalized to equivalent months of 30.4 days duration. 
The normalization of the data tends to suppress certain systematic influences 
of station elevation and latitude. 

The coastal stations of Tavoy and Rangoon, Burma shown in Figure 7a 
exhibit only a single, advectively fed surge of rainfall. All other stations 
have some continental fetch to the southwest and show evidence of an earlier 
season, evaporatively fed surge as well. The northern inland stations of 
Hanoi, Chapa and Luang Probang are the nearest to being in phase witfi tiie 
monsoon cycle exhibited at Tavoy but show some distortion attributed to early 
season evaporation from the land. Mandalay, Bangkok and Sakaerat Forest all 
show an early season, evaporatively fed rainfall maximum in May with a signifi- 
cant depression of rainfall from June to August, the months of peak rainfall 
along the coast of Burma. This suppression of rainfall is undoubtedlv due 
to the screening effect of tlie mountains and highlands at the windward reaehe 
of these stations. The monsoonal rains do not reach their maximum at these 
locations until September. Saigon and Dalat exhibit similar effects but are 
farther removed from the drying influence of the Burmese and Malayan mountain 
ranges and are also influenced by extensive flow over tlie Gulf of Thailand 


during the southwest monsoon. 


Nhatrang and Rota Bharu liotlt show significant 
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phase differences from the other locations. Kota Bharu, being on the south- 
western shore of the Gulf of Thailand with a highland area to its south and 
southwest, shows generally low rainfall during the entire period of the south- 
west monsoon but a maximum during the early months of the northeast monsoon 
when it is leeward of the gulf. Nhatrang exhibits a similar seasonal 
delay of its precipitation maximum since it has a large continental fetch 
during the southwest monsoon and is just leeward of highlands near the east 
coast of Vietnam but is subject to the maritime influence of the South China 
Sea in the later months of transition to the northeast monsoonal flow. 

Figures 8a-d show the complex monsoon cycle of Southeast Asia in a geo- 
graphic perspective. The surfaces represent the percent of annual total 
rainfall which falls in each of four particular months chosen to illustrate the 
phase relationships over the area. Kota Bharu data is not included in these 
figures due to its unique situation relative to the distribution of water 
and land. The normalization again tends to suppress some systematic eleva- 
tion and latitude effects but differences in fetch and orography contribute to 
the complexity of the distributions. Significant variations are apparent 
between the higher latitude, interior stations, which are the farthest from 
the tropical waters of the Andaman Sea and Gulf of Thailand, and the stations 
in Thailand and at Mandalay which are closer to but protected from tlie sea 
by the mountain ranges to the southwest. The opposing coastal stations of 
Tavoy, Burma and Nhatrang, Vietnam illustrate tlie dramatic phase differences 
in monsoonal rains across the region. The soutliwest monsoon seems to be best 
illustrated by the coastal stations of Tavoy and Rangoon, Burma wliere the 
monsoonal rains begin to build during May and dime reaching a maximum in 
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July. The evaporation from land areas is relatively high during the early 
months so that interior locations also show an early season surge. However, 
locations such as Mandalay, Burma and Bangkok and Sakaerat Forest in Thailand, 
being leeward of the coastal mountain ranges of Northern and Southern Burma 
and Malaysia, show a decline in rainfall during June and July, the nontlis 
of maximum rainfall on the windward side of these mountain ranges. These 
stations do not fully exliibit the monsoonal effect until August and Septenbii. 
The higher latitude inland stations of Chapa, Hanoi and Luang Prabang are 
more nearly in phase with the coastal stations of Tavoy and Rangoon. llie 
stations along the southeast coastline show a significant delay of the 
monsoonal rains with Nhatrang, for example, experiencing maximum rainfall 
in October and November when the transition from the southwesterly monsoon 
regime to the northeasterly monsoon regime impends. 

5. SUMMARY /\ND CONCLUSIONS 

The characteristics of the climate of Sakaerat Forest, Thailand, iiavo 
been discussed in tills paper. The area is influenced by a cool, dry north- 
east monsoonal flow generally off the Asian continent during the period 
from November to February and by a warm, moist southwest monsoonal I 1 ow 
off the /Indaman Sea and the Gulf of Thailand during tlie period from April 
to September. This latter period is characterized by two distinci r.iinf.ill 
surges in May and September separated by a period of lesser rainfall duriiii 
June and July. The climatic regime at Sakaerat Forest is somewhat 
to that at Bangkok, 190 km to the southwest. However, differences in <'li'v.i- 
tion and significant differences in the lieat and moisture balance in the 
forest yields a consistent cool bias of about J.5°C in the Sakaei at I', i i t 
and a more complicated interruption of tlie annual cycle of relative hui'iid.i: 
compared with Bangkok. Both stations exhibit lowering temperatures .iiul 


relative humidities during the cool northeast monsoon. The evaporation 
increases during this period due to low humidity values. The divergent 
flow around the Asian continental high pressure cell suppresses the pre- 
cipitation, however, and leads to an advective flux of moisture out of the 
Thailand region and relatively low precipitable water values. 

During the warm, moist southwest monsoon, the temperature, relative 
humidity, precipitable water and precipitation are all relatively high at 
Sakaerat Forest. Sakaerat, Bangkok and several otlier stations in southeast 
Asia show double maxima in precipitation during the season, the first 
surge during spring being supported largely by high amounts of evaporation 
but the second and larger surge during late summer and fall being supported 
by a large advectlon of moisture into the region from maritime sources 
to the south an ^ southwest. Tavoy and Rangoon, Burma, on the eastern and 
northern shores of the Andaman Sea, show single and very large monsoonal 
rainfall effects extending from May to September and peaking in July, 
about two months earlier than the surge in Thailand. Nhatrang, Vietn.un, 
and other eastern coastline stations in Southeast Asia, sliow relatively 
weak early season rainfall surges and the late season precipitation maximum 
comes one to two months later than the surge in Thailand. 

Thus, the precipitation cycles in Soutlieast Asia seem to be a complex 
function of the seasonal shift of the circulation pattern in relation to 
tlie distribution of sources of maritime air, and windward reaches of ilie 
warmed land evaporation sources. Forest cover, such as that at Sakaerat 
Forest, f\irtlier modifies the local heat, momentum and moisture budgets so 
that exchanges of these properties between the earth and atmosphere, and 


mean values of properties, are generally different than at urbanized but 
otherwise comparable locations in Southeast Asia. 


I 


6. SUPPLEMENTAL INFORMATION 

In this Section, various supplementary results are included which 
further characterize the climatological conditions at the TREND experiment 
site . 

In Section III.l, mention was made of a comparison between daily 
average values of hourly effective temperature and effective temperature 
computed from daily mean temperature and relative humidity for three parti- 
cular days. Table (8) shows the details of this comparison for December 14, 
1967; April 19, 1970; and August 13, 1970. 

Table (9) gives the effective temperature distribution for daily mean 
values. The entries are the period average values of the number of days in 
eacli month with daily mean effective temperature in the intervals shown. 

Ttie entries in parentheses are percentages of the total number of days of 
that particular month during the period of record at Sakaerat. The effec- 
tive temperature gives a reasonable first approximation to heat stress 
problems in tropical countries. It can furnish guidelines for clotliing 
to be Worn and to levels of physical labor that can be tolerated. Effective 
temperature is a fairly good measure of sultriness. This is, of course, 
a subjective sensation and, for non-indigenous individuals, a matter of 
degree of acclimatization. 
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Table 8. Daily mean of hourly effective temperatures compared with values 
determined using daily average temperature and dally average relative humidity. 


T ime 

Dec. 14 

1967 

Effective 

April 

19, 1970 

Ef f ective 

August 

13, 1970 

Effect i V 


Temp. 

R. H. 

Temp. 

Temp. 

R. H. 

Temp. 

Temp. 

R. H. 

Temp . 

0100 

13.5°C 

83Z 

13.0°C 

27. 7°C 

18% 

25. 8°C 

24.0°C 

94% 

23. 5"(: 

0200 

12.7 

92 

12.5 

27.1 

78 

25.3 

24.0 

94 

23. 5 

0300 

12.0 

86 

11.7 

27.0 

75 

2 5.1 

24.0 

94 

23.5 

0400 

11. 2 

92 

11.0 

27.0 

71 

24.9 

24.0 

92 

23.4 

OSOO 

10. 7 

88 

10.5 

27.0 

66 

24. 5 

23.0 

100 

23.0 

0000 

11.3 

86 

11.0 

26.9 

66 

24.4 

22.9 

98 

22.7 

0700 

12.0 

78 

11.6 

28.0 

63 

25.1 

25.0 

81 

13.9 

0800 

13.4 

72 

12.9 

30.0 

58 

26. 3 

26.0 

80 

24.7 

0900 

14. 7 

62 

13.9 

32.0 

56 

27.6 

28. 5 

70 

25.9 

1000 

16. 7 

56 

15.5 

34.0 

52 

28. 8 

30.5 

59 

26.7 

1100 

18. 3 

54 

16.8 

36 .*0 

46 

29.6 

31 .2 

56 

27.0 

1200 

19.5 

50 

17.5 

37.0 

38 

29.4 

33.0 

56 

2 8.4 

1300 

20. 3 

48 

18.1 

38.0 

39 

30.1 

33.0 

56 

28.4 

1400 

21.8 

47 

19.3 

36. 3 

42 

29.4 

31.9 

70 

28.9 

1300 

22.2 

47 

19.7 

36.7 

44 

29.9 

24.2 

98 

24.0 

1600 

21.2 

48 

19.0 

35.0 

48 

29.2 

25.1 

98 

24. <1 

1700 

18.0 

64 

16.7 

34.0 

48 

28.3 

25.0 

98 

24.8 

1800 

15.0 

83 

14.5 

32.1 

54 

27.5 

25.0 

93 

24.6 

1900 

14.0 

87 

13. 7 

30.2 

62 

26.8 

25.0 

86 

2 4 . ! 

2000 

13. 3 

88 

13.0 

30.3 

60 

26. 7 

25.0 

88 

24.3 

2100 

13. 2 

90 

13.0 

29.9 

66 

26.8 

24.0 

98 

23.8 

2200 

12.9 

89 

12.7 

29.2 

70 

26.6 

23.9 

100 

2 3.9 

2 300 

12. 7 

92 

12.5 

28.9 

76 

26.6 

23.2 

100 

2 3.2 

2400 

12.5 

86 

12.2 

28.2 

80 

26.6 

23.9 

98 

23.7 

Average 

15. 1 

73.7 

14. 3 

31.2 

59.8 

27.1 

26. 1 

85. 7 

24. 7 


Effective Temperature 14.6 27.4 2S.1 

determined from daily 
mean temperature and 
relative humidity 



Table <i . Effective Ter.per .ture Distribution for Dally Mean 















































Tlio data tor Sakaerat Experiment Station show tliat for persons from 
moderite latitnde couiitries tlie period from April through October is 
uaitorTiiLy in the uncomfortable" range. December and January are the only 
months when most of the days are in the "comfortable" range and some, 
proluibly not more than one out of five, are definitely cool. The area of 
fable (')) between the dotted lines represents the comfort ,";one for lightly 
clad U. S. subjects. 

on the otlier hand, the liigh stress limit of 30°ET is also rarely 
exc'et'do.d, 1 he few v'alnes that were noted were in mid-day in the pre— 
-'ummer monsoon season. In suclt conditions iieavy physical laiior iias to he 
avoided to red\ice the risk of heat stroke, 

i iib le (10) shows monthly average values of cl imato If^gic.'il elements 
avenged over the period of record at Sakaerat Experiment Station. These 
ilata are pletted as Figure A of this report. 

Data of Table 2 of this report are plotted In Figure 9. 
lubles (lla,b,c,d) give values of temperature, relative iiumidity and 
pr.'cip i tation on ,an hourly basts for the months of March, Jiuie, September 
■and Deceaniier. fliese data are plotted as Figures 5a,b,c, of this report. 
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rablelO Monthly average values of climatological elements averaged over the period of 





eriment Station. 
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11.67 4.1 0.16 

12.03 82.5 2.67 

12.47 116.4 3.90 

12.83 171.6 5.57 

13.00 122.3 4.03 

12.93 77.8 2.43 

12.63 122.2 3.90 

12.22 252.2 8.70 

11.82 120.4 3.90 

11.43 51.6 1.70 

11.25 0.0 0.00 



11.25 

12.1 


94.4 3.11 

132.5 


Evapora 


Total Dail 


124.9 3.94 

177.8 5.33 

194.2 6.27 

168.3 5.61 

145.7 4.70 

125.6 4.19 

133.0 4.29 

115.4 3.77 

92.9 3.10 

90.2 2.75 

101.2 3.38 

133.5 4.31 

133.6 4.30 
1602. 7 



- 56 . 

,8 

5 . 

.8 

156 . 

3 

00 

6 

- 53 . 

,7 

- 141 . 

2 

- 39 . 

,2 

0 

1 

2 


- 170.9 
- 76.4 
- 62.0 
1.1 
22.1 
- 56.7 
5 . 7 
157.9 
28 . 7 
- 74.5 
- 150. 3 





















TABLE 11a. DIURNAL VARIATION OF CONDITIONS AT SAKAERAT EXPERIMENT STATION 


FOR MARCH. 
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TABLE lib. DIURNAL VARIATION OF CONDITIONS AT SAKAERAT EXPERIMENT STATION 
FOR JUNE. 
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TABLE 11c. DIUBNAL VARIATION OF CONDITIONS ATSSAKAERAT EXPERIMENT STATION 
FOR SEPTEMBER. 
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In order to guide the selection of short periods for intensive study 
of energy balances, turbulent fluxes, and the like, the entire period of 
record at Sakaerat Experiment Station was scanned to pick out days and/or 
periods of significant weather conditions. 

Table (19 gives the dates of days with more than eight hours of sun- 
shine. The number of hours of sunshine is also shown as well as a coded 
symbol for the precipitation events of each day. Thus, those days with a 
blank precipitation code represent clear days with no rain during the 
twenty-four hour period; dates with (*) or (+) are clear or mostly clear 
days but with rainfall amounting to less than an inch for the twenty-four 
hour period; etc. 

Table (13) gives the dates of days with less than three hours of sun- 
shine. Again, the number of hours of sunshine and the precipitation code 
are also shown for each day. Thus, entries with blank precipitation codes 
are cloudy days but days for wtiich no precipitation fell in a twenty-four 
hour period; days with (++) are rainy days with more than an incli of pre- 
cipitation in a twenty-four hour period; etc. 

Finally, Table (14) gives dates of days with one full hour of sun- 
shine at noon during the period of record at Sakaerat Experiment Station. 



DAYS WITH MORE THAN 8 HOURS OF SUNSHINE 
DURING THE PERIOD OF RECORD AT SARAERAT EXPERIMENT STATION 
(Precipitation (P) Code: blank No rain, * between O-lOmro rain, + between 10-25inm rain, 
-H- greater than 25mm rain, M information missing) 


Date 

Hours P 

Sunshine 

Date 

Hours P 

Sunshine 

Date 

Hours 

Sunshine 

P 

Date Hours 

Sunshine 

1967 


1968 


1968 



1968 


11-17 

8.8 

1-5 

8.9 

3-5 

8.4 


5-26 

8.8 

11-18 

8.8 

1-6 

8.8 

3-6 

8.5 


5-20 

8.8 

U -19 

8.6 

1-7 

9.0 

3-7 

9.1 


6-19 

8.1 

11-20 

8.4 

1-11 

8.5 

3-10 

9.2 


6-20 

8.4 

12 -A 

8.1 

1-12 

9.0 

3-11 

8.9 


6-21 

8.4 

U -5 

8.9 

1-13 

9.1 

3-14 

9.3 


6-22 

8.9 

12-6 

8.9 

1-14 

8.5 

3-15 

8.9 


6-23 

8.9 

12-11 

9.1 

1-16 

8.9 

3-17 

9.0 


6-24 

8.3 

12 -lA 

8.2 

1-17 

8.4 

3-18 

8.8 


7-7 

8.0 

12-15 

9.0 

1-19 

8.3 

3-19 

8.4 


7-8 

8.2 

12-16 

0.9 

1-20 

8.6 

3-20 

8.6 


7-19 

9.5 

12-17 

8.6 

1-22 

9.0 

3-26 

8.7 


7-28 

8.0 

12-18 

8.1 

1-28 

8.9 

3-27 

9.3 


8-7 

8.5 

12-19 

8.9 

1-29 

9.2 

3-28 

8.5 


8-19 

9.4 

12-20 

8.7 

1-31 

9.0 

3-29 

8.3 


8-20 

9.4 

12-21 

8.8 

2-2 

9.8 

3-30 

8.7 


8-24 

8.6 

12-22 

8.7 

2-3 

9.9 

3-31 

8.0 


8-25 

9.3 

12-23 

8.9 

2-4 

10.0 

4-1 

9.7 


8-26 

8.8 

12-24 

8.7 

2-5 

10.0 

4-2 

9.8 


8-29 

10.3 

12-25 

8. 1 

2-6 

9.6 

4-3 

10.1 


9-9 

8.7 

12-26 

8.6 

2-7 

10.1 

4-4 

9.3 


9-21 

8.8 

12-27 

9.0 

2-8 

9.5 

4-15 

8.4 


9-22 

8.3 

12-28 

8.9 

2-9 

8.9 

4-20 

10.3 


9-27 

8.0 

12-29 

8.2 

2-10 

10.1 

4-27 

8.5 


9-30 

8.1 

12-30 

8 . 3 

2-11 

9.9 

5-11 

8.7 

* 

10-2 

9.7 

12-31 

8.9 

2-13 

8.9 

5-16 

8.5 


10-3 

9.1 


— 

2-15 

9.4 

5-20 

9.3 


10-23 

9.6 

1968 


2-16 

9.6 

5-21 

10.1 


10-24 

9.7 



2-18 

9.2 

5-23 

10.4 


10-25 

8.7 

1-2 

8.9 

2-27 

8.4 

5-24 

9.6 


10-26 

8.7 

1-4 

8.7 

2-29 

8.5 

5-25 

8.3 


10-28 

9.4 




Date Hours P Date Hours P Date Hours P Date Hours P 

Sunshine Sunshine Sunshine Sunshine 



10-29 

9.5 

12-15 

9.2 

2-12 

10.0 

4-21 

10.5 


10-31 

10.1 

12-19 

8.8 

2-16 

8.6 

4-23 

9.2 


11-1 

9.6 

12-20 

9.2 

2-17 

9.1 

4-24 

8.2 


11-2 

9.4 

12-21 

9.0 

2-18 

9.5 

4-27 

8.5 


11-3 

9.5 

12-22 

9.1 

2-19 

9.8 

5 - 3 

8.0 

* 

11-4 

9.5 

12-23 

9.2 

2-20 

9.7 

5-5 

9.1 


11-5 

9.4 

12-24 

9.0 

2-11 

9.8 

5-9 

9.8 


11-7 

9.3 

12-25 

9.0 

2-22 

9.7 

5-12 

8.4 


11-8 

8.5 

12-26 

8.4 

2-23 

9.3 

5-20 

8.1 


11-11 

9.5 

12-27 

8.3 

2-24 

8.0 

5-21 

8 . i 


11-12 

9.3 

12-28 

8.5 

2-25 

8.0 

6-2 

8.4 


11-13 

9.4 

12-30 

8.9 

2-28 

8.0 

6-29 

9.0 

* 

11-16 

9.0 

— 


3-2 

8.7 

6-30 

8.0 

* 

11-17 

9.0 

1969 


3-3 

8.6 

7-5 

9.5 


11-18 

9.3 



3-4 

8.9 

7-14 

9.1 

* 

11-19 

9.0 

l-I 

8.6 

3-5 

8.7 

7-16 

9.5 

* 

11-20 

9.1 

1-2 

9,1 

3-8 

9.3 

8-14 

9.5 

* 

11-21 

9.0 

1-3 

9.2 

3-10 

9.5 

8-16 

9.0 


11-22 

9.2 

1-8 

8,4 

3-20 

8.2 

8 - 19 

8.8 

* 

U -23 

9.0 

I-IO 

8.8 

3-27 

9.4 

8-27 

9.4 

* 

11-24 

8.1 

l-ll 

8.4 

3-28 

8.5 

8 - 31 

9.0 


11-26 

8.0 

1-20 

8.4 

3-29 

8.5 

9-8 

8,6 


11-28 

9.0 

1-21 

8.8 

3 - 30 

8.7 

10-11 

8.5 


12-2 

9.1 

1-23 

8 . 1 * 

4 - 3 

H.h 

10-16 

8.0 


12-3 

8.9 

1-26 

8.2 

4-6 

9.0 

10-17 

8.5 

* 

12-4 

9 . 1 

1-27 

9.5 

4-7 

8.9 

10-18 

9.0 


12-5 

9.1 

1-29 

9.3 

4-8 

9.5 

10-19 

H . 1 


12-7 

9.0 

2-5 

9.0 

4-9 

8.7 

10- :o 

8.0 

* 

12-8 

9.2 

2-6 

8.6 

4-1 3 

9 . 1 

10-28 

H.U 


12-9 

8.6 

2-7 

9.5 

4-14 

9.4 

10-29 

H .5 


1 11 

9.1 

2-8 

9.7 

4-15 

10. 0 

10 - 30 

8.0 

■4 

12-12 

9 . 1 

2-9 

9.8 

4-16 

8.6 

1 1-9 

H .8 


I 2 -n 

8.9 

2-U) 

9,4 

4-17 

9.7 

ll-lO 

8 . 3 


tJ- 14 

9.0 

2-11 

9.7 

4-18 

9.4 

n-ii 

‘^6 



■ 


r ^ 

1 

i 

I 

TABLE 12 (CONTINUED) 


Date 

Hours 

Sunshine 

P 

Date 

Hours P 

Sunshine 

Date 

Hours 

Sunshine 

P 

Date Hours 
Sunshine 

P 

1969 



1970 


197^ 



1970 




11-12 

3.4 


1-3 

10.0 

2-24 

9.7 


4- 30 

10.0 


11-13 

8.7 


1-7 

8.5 

2-25 

9.7 


5-2 

8.9 


U-20 

8.3 


1-8 

9.0 

2-26 

9.8 


5-3 

8.3 

A 

11-21 

8.8 


1-9 

9.1 

2-27 

10,2 


5-8 

9.4 


11-22 

8.7 


1-10 

9.1 

2-28 

10.1 


5-26 

9.3 


11-27 

9.1 


1-11 

8.6 

3-2 

8.3 


5-27 

9.1 


11-30 

9.1 


1-12 

10.0 

3-3 

8.4 


6-2 

8.7 


12-1 

9.0 

M 

1-13 

9.0 

3-4 

8.8 


6-10 

9.4 


12-2 

8.7 

M 

1-14 

8.4 

3-5 

8.2 


6-16 

8.5 


12-3 

8.9 

M 

1-18 

9.0 

3-9 

8.6 


6-24 

9.2 

A 

12-4 

8.9 

M 

1-19 

9.3 

3-10 

8.1 


6-29 

8.3 

A 

12-5 

8.9 

M 

1-20 

8.1 

3-11 

8.9 


7-2 

8.0 


12-6 

8.7 

M 

1-21 

8.4 

3-12 

B.8 


7-3 

9 . 6 


12-7 

8.9 

M 

1-23 

9.2 

3-13 

8.9 


8-9 

9.1 


12-8 

8.9 

M 

1-24 

8,7 

3-14 

8.9 


8-10 

9.4 


12-10 

9.1 

M 

1-26 

8.7 

3-15 

9.1 


8-30 

10.2 


12-11 

9.3 

M 

1-27 

9.3 

3-16 

9.1 


8-31 

8.3 


12-12 

9.2 

M 

1-29 

9.3 

3-17 

8.8 


9-4 

8.5 

A 

12-14 

8.6 

M 

1-30 

8.9 

3-18 

8.8 


9-8 

9.2 


12-18 

8.1 

M 

1-31 

9.5 

3-25 

8.8 


9-25 

8.3 

A 

12-21 

8.9 

M 

2-1 

8.7 

3-28 

8.3 


10-4 

9.1 

M 

12-22 

8.9 

M 

2-2 

9.5 

3-29 

8.8 


10-5 

8.4 

M 

12-23 

8.9 

M 

2-3 

9.4 

3-30 

8.6 


10-8 

8.6 

M 

12-24 

8.9 

M 

2-4 

8.5 

4-9 

8.7 


10-10 

8. 1 

M 

12-26 

9.2 

M 

2-5 

9.5 

4-10 

9.6 


10-20 

8.5 

M 

12-27 

8.8 

M 

2-6 

9.5 

4-14 

9.6 


10-25 

8.5 

M 

12-29 

9,2 

M 

2-7 

10.0 

4-15 

10.3 


11-5 

9.3 

M 

12-30 

9.9 

M 

2-8 

8.9 

4-16 

10.5 


11-6 

9.0 

M 

12-31 

8.8 

M 

2-9 

9.7 

4-17 

9.9 


11-7 

9.2 

M 




2-19 

9.6 

4-18 

10.6 


11-12 

9.3 

M 

1970 



2-20 

10.7 

4-19 

10.1 


11-13 

8.8 

M 




2-21 

10.2 

4-20 

10.4 


11-14 

8.3 

M 

1-1 

9.0 


2-22 

9.9 

4-22 

8.1 

* 

11-15 

8.4 

M 

1-2 

9.0 


2-23 

9.8 

4-29 

8.0 

A 

11-17 

8.5 

M 




Date 

Hours 

Sunshine 

P 

1970 

11-20 

8.3 

M 

11-21 

8.1 

M 

12-11 

8.0 

M 

12-16 

8.8 

M 

12-17 

8.8 

M 

12-18 

8.5 

M 

12-19 

8.3 

M 

12-20 

8.6 

M 

12-21 

8.7 

M 

12-25 

8.8 

M 




I 


table Li. days with less than 3 HOURS OF SUNSHINE DURING THE PERIOD OF RECORD AT SAKAERAT EXPERIMENT STATION 
(Precipitation (P) Code: blank, no rain, * between O’-lOimn, between lO-25nun, -h- greater than 25inin, M information missing) 


Date 

Hours 

Sunshine 

P 

i^b/ 

11-27 

0.9 


11-28 

2.7 


12-8 

1.8 


1968 

2-2A 

2,7 

* 

3,2 

0,1 


4-S 

2,3 


4-19 

1.7 

* 

4-26 

2.2 


4-30 

2.8 

■H* 

3-1 

0.9 

A 

3-4 

2.5 

A 

3-6 

2.7 

+ 

3-8 

2.8 

A 

3-29 

0.1 


6— 1 

1.2 

•f 

6-2 

0.0 


6-3 

U.l 

A 

6-4 

1.7 

+ 

6-5 

0.2 


6-8 

1.7 

A 

6-17 

1.1 


6-2 7 

1.6 

A 

6-29 

1.1 

A 

6-30 

1.3 

A 

7-5 

1.8 


7-22 

2.4 

A 

7-23 

1.2 


7-26 

0.9 


7-27 

2.2 

A 


Date 

Hours 

Sunshine 

P 

1968 

8-3 

0.9 


8-4 

2.1 

A 

8-5 

0.4 

+ 

8-10 

0.1 


8-U 

0.0 


8-12 

0.4 

A 

8-13 

0.0 


8-14 

0.6 

A 

8-18 

0.2 

A 

8-21 

2.3 


9-2 

2.8 


9-5 

0.4 

A 

9-6 

0.0 


9-7 

0.0 


9-12 

2.4 

-H- 

9-13 

0.3 

+ 

9-28 

2.7 

A 

10-10 

2.2 

+ 

10-15 

2.9 

A 

10-16 

2.0 

A 

10-21 

0.3 

A 

10-22 

0.1 

++ 

11-25 

0.5 


12-6 

1.9 


12-10 

1.8 


1969 

1-6 

0.0 

A 

1-25 

1.9 

A 

2-2 

1.6 


3-14 

1.7 





Date 

Hours 

Sunshine 

P 

Date 

Hours 

Sunshine 

P 

1970 



1970 



1-4 

1.1 


8-3 

0,0 

A 

1-5 

0.1 


8-14 

0.4 


1-6 

l.l 


8-18 

0.0 

A 

1-15 

2.5 

+ 

8-19 

2.7 


1-16 

1.2 


8-20 

0.0 

A 

1-17 

2.4 


8-23 

0.0 

A 

2-11 

1.3 


8-25 

1.1 

■f 

2-12 

0,0 

A 

8-26 

0.0 

A 

3-6 

2.7 


9-2 

l.l 


3-23 

0.0 


9-9 

1.5 

A 

3-26 

2.3 


9-10 

1,1 

+ 

4-4 

2.6 

A 

9-11 

2.1 

A 

4-13 

1.5 

M- 

9-17 

1.3 

++ 

4—26 

1.4 

A 

9-27 

2.8 

M 

3-15 

2.6 

-*-*■ 

9-29 

2.8 

M 

3-16 

1.5 


10-11 

1.5 

M 

5-23 

1.4 

*■ 

10-12 

2.9 

M 

6— H 

1.9 

■M- 

10-16 

0.0 

M 

6-13 

2.2 


10-19 

0.4 

M 

6-14 

1.0 

A 

U)-26 

0.1 

M 

6-15 

2.4 


10-27 

2.1 

M 

6-21 

0.5 

A 

10-29 

1.7 

M 

6-27 

2.5 

A 

10- 10 

0.0 

M 

6-28 

1.3 

A 

10-31 

0.0 

M 

7-5 

1.7 

A 

11-1 

0.6 

M 

7-H 

2.5 


ll-H 

0.8 

M 

7-9 

1.4 

A 

ll-io 

1 .6 

M 

7-10 

2.1 

A 

ll-ll 

0.2 

M 

7-14 

0,4 

A 

11-29 

0.0 

M 

7-15 

0.0 


11-30 

0.7 

M 

7-16 

0.4 


12-1 

0.5 

M 

7-17 

0.0 

A 

12-1 

2.1 

M 

7-18 

0.2 

A 

12-4 

2.9 

M 




Date 

Hours 

Sunshine 

P 

Date Hours 
Sunshine 

P 

1969 



1969 



3-17 

1.8 

-H- 

8-7 

0.0 

A 

3-18 

0.0 

A 

8-8 

0.0 

A 

4-5 

0.0 

++ 

8-9 

0.2 

A 

4-26 

0.0 

+ 

8-10 

0.0 


4-28 

1.8 

+ 

8-11 

0.7 


5-7 

0.0 


8-12 

0.0 

A 

5-16 

2.6 

-H- 

8-11 

2.0 

A 

5-27 

1.4 

A 

9-3 

0,0 

A 

5-29 

0.0 

A 

9-4 

0.0 

+-f 

5-30 

0.2 

A 

9-5 

2.8 

A 

5-31 

0.8 

A 

9-6 

1.8 

+ 

6-5 

1.3 

A 

9-15 

2.1 

A 

6-3 

1.5 


9-16 

0,1 

+A 

6-12 

2.3 

A 

9-17 

0.5 

A 

6-13 

1.7 

+ 

9-20 

1.0 

-4-6 

6-18 

0.0 


9-21 

2.7 

+ 

6-25 

0.6 

A 

10-2 

0.0 

-M- 

7-3 

0.8 

A 

10-5 

1.8 

A 

7-8 

0.0 

+ 

10-6 

1.5 


7-11 

1.5 

A 

10-7 

0.6 


7-12 

0.0 

A 

10-9 

2.5 

A 

7-13 

0.2 


10-27 

1.6 

A 

7-18 

2.4 

+ 

11-1 

1.2 

A 

7-19 

0.0 

+ 

11-2 

0.7 

A 

7-20 

1.5 


11-3 

1.3 

AA 

7-23 

0.3 

A 

11-25 

2.5 


7-24 

2.7 


11-29 

2.9 


7-26 

2.8 


12-19 

1.6 

M 

7-28 

0.0 

A 

12-28 

1.5 

M 

7-29 

0.0 

A 




7- 30 

0.5 





8-5 

0.7 






Date 

Hours 

Sunshine 

P 

1970 

12-5 

0.0 

M 

12-10 

2.4 

M 

12-14 

0.1 

M 

12-22 

2.5 

M 

12-29 

2.8 

M 


I 


TABLE 14. Days with One Full Hour of Sunshine at \oon 
during the Period of Record at 
Sakaerat Experiment Station 
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SECTION IV - ANALYSIS OF V/IND AND TEMPERATURE 
PROFILE CH;\RACTERISTICS 

1. PRELIMINARY ,VNA1.YSIS OF TEMPERATURE AND WIND PROFILES 

Temperature and wind speed data from the two towers were subjected 
to preliminary analysis for tliree different periods of the year. Figures 
3a,b,c show typical temperature and wind speed profiles for each tower 
during day and night periods for January 1-lA, 1970, during the cool, 
dry northeast monsoonal flow; for June 20-30, 1970 and September 6-12, 

1970 during the warm, moist southwest monsoonal flow. The profiles repre- 
sent averages of all available profiles during the hours shown and over 
the periods mentioned above. 

The cleared area generally maintains higher temperature and wind 
speed than the forest. There is an inversion of temperature just beneath 

•k 

the forest canopy in January during the daytime with strong lapse above. 

The warmest temperature is at 30 m, generally near the canopy top. This 
warm layer tends to carry over into the clearing. Even though the tempera- 
ture in January is lower than in June or September, this warming effect 
is more pronounced in January. This is because the cloud cover is less 
during the cool northeast monsoon than during the moist southwest monsoon 
and the sun shines on the forest for longer periods of the day in January 
tlian in June or September. [See Section III of this report.] The noc- 
turnal profile of temperature for January also shows prominent radiational 
cooling at the canopy top level due, in part, to the low values of prcci- 
pitable water in the air column above the forest region during the northeast 
monsoon . 

* Baynton, et al. (1965) found similar behavior above a tropical forest. 





a) JANUARY b) JUNE c) SEPTEMBER 
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The temperature of the forest floor in all three months is consider- 
ably cooler during the day than the surface in tiie clearing. 

During the January daytime period shown, the surface temperature in 
the clearing was cooler than at 1 meter above tlie surface due to low 
shrub growth whicti was subsequently cleared — the daytime temperature 
profiles in the clearing for June and September show maximum temperatures 
at the Surface. Consistent with this, the surface temperature in the 
clearing during the night in January remained warmer than at 1 meter. 

This is due, in part, to the displacement of the effective level of 
radiational cooling to the top of tlie shrub growth and possibly also to 
a drainage cooling effect above the surface of the sloping terrain which 
is evident in the profiles for both the forest and the clearing. During 
June and September, tlie higlier moisture content of the air tends to 
suppress any strong surface radiational cooling in the clearing during 
the nighttime hours. This leads to a diurnal change in the temperature 
profile which amounts to more of a shift to cooler temperatures at all 
levels during the night than a lapse reversal at low levels. 

Tlie nocturnal temperature profile in the clearing during the drier 
month of January also shows a rather strong inversion above 1 meter 
while the strongest Inversions in tlie forest occur above the lov/er and 
upper stories of the forest. The nocturnal profiles of temperature 
in the clearing during June and September do not exhibit the prominent 
inverted structure of January. 



The day-night temperature difference at 30 m, near tlie top of the 
forest canopy, is 4°C for all three months while at the surface the 
difference is only about 1°C or less. In contrast to this, in the clearing 
the differences are generally largest at the surface. Of course, such 
differences between six-hour averaged temperatures under-estimates the 
true diurnal variations in temperature. For example, In Section III of 
tills report, we have established the mean diurnal variation of air 
temperature at Sakaerat Forest to be about 7°C in June, 8-9°C in 
September and 11-12°C in December. 

The sub-surface temperature profiles show the clearing to remain 
always warmer than the forest soil even to a depth of 1 meter. The 
surface temperature undergoes a stronger diurnal variation in the clearing 
as expected. The temperature differences between the floors of the 
forest and clearing during the daytime period are about 3°C, 3.5°C and 5.3°C 
in January, June and September respectively while at night the diffe- 
rences are only about l^C or less for all three months. At one meter 
depth, the differences in temperature between tlie forest and cleared 
area are about 1°C, 2.5°C and 2°C for January, June and September 
respectively with very little variation from day to night. 

The wind speed profiles show the drag effect of the forest canopy 
in all seasons and at all times of the day. Above the forest canopy, 
the wind speed tends to show a logarithmic behavior while within the 
canopy, the wind speeds are low with rather unremarkable vertical 
variation. The anemometer at the 40 m level of tlie tovier in the forest 
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had a bearing failure during January so that the values shown for 40 m 

* 

are not representative of the winds at that level. This anemometer v:as 
replaced shortly after the two week period represented by the January 
data and does not contaminate the profiles for June or September. 

The day-night differences in the wind speeds indicate a nocturnal 
Increase of wind speed at all levels in June and September during the 
southwest monsoon, but a daytime maximum of winds in January during the 
northeast monsoon. This is the case in both the forest and clearing. 

The wind speed profiles in the clearing depart considerably from a 
logaritlimic relation indicating that tlie cleared area is significantly 
influenced by the . urrounding forest. It should be noted, however, 
tliat Che data appear to define two distinguishable layers with the wind 
speed in Che first 10 or so meters behaving differently than the data 
above. Tills is particularly evident in the June and September periods. 
Presumably, the wind cliaracterist ics at low levels are closely influenced 
by roughness and stability conditions near the floor of tlie forest wliile 
tlie forest canopy influences the wind beliavior at upper levels along 
the tower. 


* Personal communication with J. Zabransky, Field Experimenter during tlie 
TREND project. 
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2. S tability Characteristics of the Forest Area 

For an assessment of the mean stability characteristics of the 
Sakaerat Forest area, values of the Richardson number were computed for 
several layers using wind and temperature profiles averaged for each 
thirty minute period of the day and further averaged over the periods in 
January, June and September represented by Figure 10, Figure 11 shows 
the diurnal variation of the various Richardson numbers for these three 
months. Richardson nuiabers are generally positive during nighttime 
hours clianging to negative values during the day. The transition times 
are about 0700 and 1700. The major exception to this cycle is during the 
month of January where the daytime warm layer development above the canopy 
mentioned previously yields unstable conditions in the layer 32-36 m, but 
stable daytime conditions in the layer 36-46 m. The two layers show more 
consistency during Jirn md September. The reader should take care to 
note the different scales in Figure lla and Figures ilb,c. 

Bulk Richardson gradient numbers v;ere computed in the manner of 
Lettau (1957) for several conglomerate layers in the forest and clearing. 
The bulk Richardson gradient number is computed from local Richardson 
numbers as follows: 

(Ri)|^ = Z (Ri)j / Z Z. 




FIGURE 11. niURNAL VARIATION OF 
RICHARDSON NUMBERS FOR VARIOUS 
LAYERS ABOVE THE FOREST c:AN01’Y 
DURING a) .lANUARY b) JIIHF c) SEl’TEMBKK 



r 
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where (Ri)j is the local Richardson number at a given geometric mean 
height Zj and the sum if over the various sub-layers comprising a 
thicker layer whose stability is to be characterized. Lettau presents 


the stability classification shown in Table 15 for layers near the surface. 
Figure 12 shows values of (ki)j^ for the layer 32m - 46 m above the forest 
canopy. The summation is carried out over the three geometric mean 
heights of 33.9 m, 40.7 m and 38,4 m. All values shown fall within the 
neutral classification of convective stability according to Lettau' s 
classification. The values for January are all positive and exhibit 
an irregular diurnal variation due to the more complicated thennal 
structure as exemplified in Figures I'Ja and 11a. The values for June 
and September are smaller in magnitude and quite consistent with slight 
positive values at night and negative values during the day from about 
0700 to 1700 local time. 

Figures 13a, b show bulk Richardson gradient nimibers for the laver 
1-16 m in tlie clearing. Geometric mean heights for tliese computations 
are 1.41 m, 2 m, 4 m and 8 m. Note that the scale for January is two 
orders of magnitude larger than that for June and September corresponding 
to much greater extremes of stability and instability during the clear 
sky conditions of January. Tlie diurnal variation during January is also 
much more irregular than during June or September although all three 
months show negative values during tiie day and positive values at night. 


Figure iJc shows bulk Ricliardson gradient numbers for the layer 16-4b m 
in tile clearing [geometric mean heigiits for these computations are 13.9 m. 



TABLE 15. STABILITY CLASSIFICATION OF BULK RIQIARDSON GRADIENT NUMBER 


(After Lettau, [1957]) 


Stability 

Class 

(Ri)^,(10 ^/m) 


Stability 

Class 

(Ri)j^,(10 

L3 =Extreme Lapse 

< than -15 

11 

=Weak Inversion 

5 to 6 

L2 =Moderate lapse 

- 15 to -10 

12 

=Moderate Inversion 

7 to 9 

LI = Weak, lapse 

9 to - 5 

13 

=Strong Inversion 

10 to 19 

N = Neutral 

- A to A 

lA 

=Extreme Inversion 

> 19 



12. Ijiill; II LcIi.'.ril'On GraditMiL Nui'bcM' fur the Microlayur 
Al)ov(' till' )''orL;.‘;f (l.inupy fur vlniiuary, .lunc ar.d Supt unl'ur 






i 



73 


22.6 m, 33.2 m, 38.4 m]. The scale is the same as Figurel3b but, again, 
two orders of magnitude smaller than Figure 13a. Thus, the range of 
stability is much greater in the lower layer than the upper layer during 
January but comparable during June and September. The bulk, Richardson 
gradient numbers are in the neutral category of Lettau from midnight to 
about 0730 during June and September but in an inversion class during 
these hours in January. January and June siiow a shift to negative values 
during the day but September exhibits stable conditions during the daytime 
hours. All months show neutral or stable values after about 1700 local 
t ime . 

3. Roughness Characteristics of the Forest Area 

Estimates of roughness parameter, z^, datum displacement height, d, and 
friction velocity, u^, were made for the forest canopy using a large 
sample of 30-minute averaged wind profiles occurring during near-neutral 
conditions in January, June and September. Computation of (ki)j^ were 
first made for each individual 30-minute averaged wind and temperature pro- 
file for the periods January 8-lJ, 1970; June 20-26, 1970; and September 9-11, 
1970 for the layer 30-46 n just above tlie forest canopy. Tlie resulting 
values were used to screen the sample to extract only those profiles falling 
into the "N" category of Table 15. Since there was doubt about the wind speed 
measurements at 40 m in January, tlie data at tliis level were excluded for 
this period. For comparison, computations of (Ri)| were also made for the 
layers 1-16 m and 16-46 m in the clearing for the June and September periods. 



’=■■■ - ■ - - 



A frequency distribution of the number of occurrences in each stability 
class is given in Table 16. Nearly all the cases sampled for the layer above 
the forest canopy during the moist, southwest monsoon period (i.e. those 
for June and September ) fell into the near-neutral classification of Lettau. 
For the same layer in January, during the cool, dry northeast monsoon, about 
41% of the cases exhibited near-neutral conditions with slightly less than 
10% falling in various of the unstable categories and slightly less than 
half falling in various of the stable categories. For the clearing, the 
lower layer exhibits a wider variation of stability than does the upper 
layer for both June and September. Sliglitly less than 80% of the cases in 
June and about 69% of the cases in September show near-neutral conditions 
in the upper layer of the clearing; evidence of the modifying influence of 
the surrounding forest as noted earlier. 

For near-neutral conditions, it is hypothesized that tlie wind speed 
variation with height above the forest canopy is represented by the logarithmic 
relat ion 

u(z) = (u./k) In [(z - d + z ) /z ] 

" o o 

where u(z) is the wind speed, u^ is the friction velocity appropriate to 
the canopy, k * Von Karman's constant, z^ is the canopy roughness parameter 
and d is the datum displacement height. Of tlie various methods for esti- 
mating z^, d and u^ which appear in the literature, the method adopted in 
tills study is that described by Steam's (1970), although comparisons will 

be made with other methods. in Steam's method, the parameters z . d, u, 

o’* 



TABLE 16 


FREQUENCY DISTRIBUTION OF CONVECTIVE STABIl.ITY CLASSES 
FOR PERIODS IN JANUARY, JUNE AND SEPTEMBER 


Period 

No. of 


Layer 

Cases 

IT 

Jan. 8- 10, 19 70 

Fores t- (30- 46m) 

144 

3 

June 20-26,1970 
Forest- (30- 46m) 

336 

0 

June 20-22,1970 

Clearing- (l-16m) 

138 

31 

June 20-22,1970 

Clearing- (16-46m) 

138 

1 

Sept. 9-11, 1970 
Forest- (30-46m) 

144 

0 

Sept. 9-11, 1970 

Clearing- (l-16m) 

144 

38 

Sept. 9-11, 1970 

Clearing- (16-46m) 

144 

0 


Stability Class 


hi 

U 

N 


12 

U 


2 

9 

59 

6 

6 

6 

52 

0 

0 

336 

0 

0 

0 

0 

8 

12 

71 

15 

0 

0 

1 

7 

16 

108 

3 

0 

3 

0 

0 

1 

140 

3 

0 

0 

0 

4 

7 

17 

21 

14 

15 

28 

0 

0 

99 

6 

8 

5 

26 
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are determined by an iterative procedure designed to minimize the mean 

2 

square deviation Z where 


= u(Zj^) - (u^/k) In [(z^ - d + z^)/z^] 

The iteration procedure of Steam's was applied to the ensemble of 30- 
minute averaged wind speed profiles in the layer 30-46 ra above the forest 
canopy which fell in the "N" category of Table 16. The June data were 
furtlier partitioned into a batch (A) of 201 cases for which the solar 
radiation impinging on the forest canopy was relatively steady, and a batch 
(B) of 72 cases during periods of comparaiively rapid fluctuations of solar 
radiation. The iterative procedure converged for 273 of the 336 June 
cases but only 40 of 59 cases in January and 67 of 96 cases during September. 
The cases for which convergence was not obtained generally corresponded to 
values of (Ri)^ which were negative and somewhat near the cut-off of the "H" 
classification of Table 15 suggesting that the limits of this cateogry 
should be narrower than as given in Table 15. Other investigators [for 
example, henschow and Johnson (1968)] have chosen to narrow the limits of 
tl\e near-neutral category. 

Tlie resulting estimates of tfie averages and standard deviations of the 
parameters z^, d, u^ are presented in Table 17. Tlie roughness parameter 
varied considerably over the periods examined. The mean value for the Jamiarv 
days was only 0.83 meters witli a standard deviation more than lialf of tlie mean 
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TABLE 17. AVERAGE VALUES AtlU STANDARD DEVIATIONS OF ROUGHNESS 

PARAMETER, DATUM DISPLAGEMENT HEIGHT, FRIGTION VELOCITY 
AND CANOPY STRESS FOR THE SAKAERAT FOREST. 


Period 

No. of Convergent 
Cases 

Zo(tn) 

[o^ (m) ] 

0 

d(ra) 
[aj(m) ] 

u^(m-s ^) 

[a (m-s~^)] 

t(u^) (dyn 

January 8-10,1970 

0.83 

29.53 

0.33 

1.31 

40 cases 

[0.47] 

[ 0.05] 

[0.09] 


June 20-26, 1970(A) 

5.48 

27.18 

1.42 

24.2 

201 cases 

[1.66] 

[ 0.69] 

[0.45] 


June 20-26,1970(3) 

5.62 

27.16 

1.26 

19.1 

72 cases 

[1.77] 

[ 0.72] 

[0.45] 


Sept. 9-10, 1970 

3.61 

28.09 

1.07 

13.7 

67 cases 

[1.10] 

[ 0.75] 

[0.48] 


Average (weighted 

4.69 

27.58 

1.21 

18.97 


by No. of 
cases) 
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value. The mean profiles of Figure ID indicate lower wind speeds during 
January than during June and September. For June, the mean value of forest 
roughness parameter is about 5.5 m with standard deviation only a third of 
the mean value. The values of datum displacement height show more consistency 
over the three periods but are larger in January than in September or June. 

The friction velocity is smallest in January and largest in June. 

Our results using Steam's method were compared against various empirical 
formulae which have appeared in the literature. Stanhill (1969) gives a 
relationship between d and canopy height, h, as follows: 


log(d) = (O.9793)log(h) - 0.1536 


This relationship seems to be in excellent agreement with values reported 
in the literature for vegetation canopies of lower height. Szeicz, et al. 
(1969) found an empirical relationship between z^ and h as follows: 


log(z^) = log(h) - 0.98 

Figure lA shows estimates of d and z^ using these formulae for canopy heigiiL 
vaiies of 29 m - 35 m along with our results using Steam's method and a 
siumnary of several other studies given by Leonard and Federer (1973). 

The Stanhill formula yields a displacement height some A - 10 m smaller than 
the values estimated using Steam's method. The Szeicz formula yields a 
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result which is consistent with the profile estimate for September but not 
for January or June. 

A test was made for possible correlations between z and d with wind 

o 

speed at the top of the tower (46 m) . Rauner (1960) found that z in- 
creases and d decreases with increasing wind speed while Tajchmann (1967) 
found that z^ decreases and d remained constant as wind speed increased. 

Alle-' (1968) and Belt (1969) found no significant dependence on wind speed 

if d were allowed to be different for each computation of (z , d) but when 

o 

a fixed value of d was used for an ensemble of profiles, z was found to 

o 

decrease as wind speed increased but only for a particular wind direction. 

In the present study, the profiles were not partitioned according to v/ind 
direction, however, d and z^ were computed independently for each profile. 

Table 18 shows the results of a simple correlation-regression analysis 
applied to the ensemble of 384 cases in January, June and September, It 
should be noted that the range of wind speeds above the forest canopy is not 
very large. Nonetheless, the roughness parameter xS found to be an increasing 
function of wind speed and the datum displacement height a decreasing function 
of wind speed, which supports the findings of Rauner. Both parameters 
correlate best with the wind speed near the canopy (i.e. u(32)) but the 
correlation between d and u(46) is also significant at the 5% level. 

The correlation between roughness and wind speed is sligiitly better wiien t e 
logarithm of z^ is used. 
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TABLE 18. CORRELATION AND REGRESSION STATISTICS FOR 384 CASES OF 30-MINUTE AVERAGED 
WIND PROFILES IN JANUARY, JUNE AND SEPTEMBER ABOVE THE SAKAERAT FOREST. 



5% level 



Finally, a preliminary analysis was made of a small sample of wind 
speed profiles in the clearing to determine to what extent these profiles 
can be fitted by simple algebraic relationships. Profiles of 30-minute 
averaged wind speeds for four days in January and two days in June were 
fitted, in the least squares sense, to linear and quadratic functions of 
height. Table 19 shows the results of these trials. The data at the top 
6 or 7 levels in the clearing seem to be adequately fit by a linear relation 
and the winds at all levels seem to be adequately fit by a quadratic 
relation. 

4. Summary 

An analysis of the characteristics of wind and temperature distributions 
in a tropical evergreen forest in Thailand has been presented in this paper. 
Temperatures are generally cooler and wind speeds lower during the northeast 
monsoon season as compared with the southwest monsoon period. Conditions are 
more variable during January since the skies are clearer, the moisture content 
of the atmosphere above the forest is less and the radiational budget is more 
dynamic. The atmosphere just above the canopy therefore exhibits a much 
more complex thermal character during January than June or September with a 
large diurnal variation of stability which is also highly dependent on the 
proximity to the canopy. 

Conditions during the warm, moist southwest monsoon generally show a 


much weaker diurnal variation. Cloudiness and atmospheric moisture are 
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TABLE 19. PRELIMINARY RESULTS OF FITTING WIND SPEED DATA IN 
THE CLEARED AREA TO LINEAR AND QUADRATIC PROFILE 
EQUATIONS FOR JAJJUARY AND JUNE 


Date 

Number of tower 
levels used starting 
from the top 

Number of 
Profiles 

Degree of 
Polynomial 

Explained 

Variance 

June 

20-21,1970 


6 

96 

1 

98.7 

June 

20-21,1970 


7 

96 

1 

98.1 

June 

20-21,1970 


8 

96 

1 

96.7 

June 

20-21,1970 


8 

96 

2 

98.2 

Jan. 

4-5,1970 


6 

96 

1 

79.5 

Jan. 

4-5,1970 


7 

96 

1 

79.8 

Jan. 

4-5,1970 


6 

96 

2 

82.9 

Jan. 

4-5,1970 


8 

96 

2 

83.9 

Jan. 

6-7,1970 


6 

96 

1 

97.9 

Jan. 

6-7,1970 


7 

96 

1 

97.5 
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greater during this season which tends to dampen the diurnal cycle in the 
radiation budget and, hence, the variation of stability above the canopy. 
Tliere is a much higher frequency of near-neutral conditions above the canopy 
during this season than during the cool, dry northeast monsoon season. 

For near neutral conditions in both monsoon periods, the winds in the 
microlayer above the forest canopy show a reasonable fit to a displaced 
logarithmic profile. The datum displacement height has an average value of 
about 27.6 meters but varies over a range of about 2(1/4) meters from the 
northeast to southwest monsoon periods. This seasonal variation is about 3 
times larger tlian the standard deviation of the short period sample estimate 
during the southwest monsoon period and nearly 50 times larger than that 
for the northeast monsoon sample. 

The roughness parameter for the canopy exhibits a strong seasonal 
variation with a mean value of 0.83 m for the January sample, 5.55 m for the 
June sample and 3.61 m for the September sample. The seasonal variation 
is about 2.7 times larger than the standard deviation in during the short 

period samples of June but about 10 times larger than that for January. 

_2 

The canopy stress has a mean value of 1.3 dynes-cm for January and tbie 

-2 -2 

much larger mean values of 21.6 dynes-cm for June and 13.7 dynes-cm for 
September. 

It is interesting to note that for the January cases for whicli the 
stability of the air above the canopy is much more variable and tlie wind 
speeds are smaller, the datum displacement height has its maximum value and 
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greatest consistency while the roughness parameter has its lowest value and 
least (relative) consistency. A pooled sample of the 30-minute averaged 
profiles of wind speed occurring in near-neutral conditions in January, 

June and September (some 384 cases) reveal a significant positive correla- 
tion between roughness parameter and wind speed and a significant negative 
correlation between datum displacement height and wind speed. 

5. Supplemental Information 

A. Thermal Characteristics for January 

In discussing the stability of the air layer above the forest canopy 
during the January period in Part 2 of this Section, it was noted that the 
daily range of Richardson numbers was quite large. This was attributed to 
a more dynamic variability of thermal stability. To substantiate this, the 
thermal characteristics of the layer alone may be examined. 

Figures (15) and (16) show time-height sections of temperature along 
both towers averaged over each 30-minute period of the day and over the 
days 1-14, January, 1970. The effect of the forest is evident in Figure 
( 15 ) where an inversion just below the tree tops and lapse condition just 
above are established and persist throughout the morning and early after- 
noon. No such feature is evident in the clearing area. A weak early 
morning inversion in the trunk space of the forest tends to break down 
around 1000 local time but becomes reestablished in stronger fashion during 
the afternoon. At similar levels in the clearing, lapse conditions pre- 
vail during this afternoon period. Similar features are apparent in 
Figure (17) for 19-30 June, 1970 except that the perturbation at the level 
of the tree tops seemed to be weaker during the early morning period. 
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Figure (18) shows the static stability a= [(g/c^) - (AT/.'.z)]/T computed 

for various layers defined by the top 6 levels of the forest tower. The 

values for the individual layers are most consistent during the nighttime 

hours 1730-0730 local tine. During this period the values o(25-30), ^(36-40), 

-4 -1 

a(40-46) are all within about 1x10 ra vjith j(25-30) consistently positive 
and a(36-40) and j(40-46) consistently negative. 

During the day the stabilities undergo a rather well defined variation. 
For example, consistent with the warm layer development above the canopy 
during the daytime, j ( 25-30) is positive in the stable lower portion of the 
layer with j(30-32) and a(32-3b) highly negative in the unstable upper portion 
of tliis layer. In the higher layers near the top of the tower, '(40-46) and 
j(36-40) both show large stable values during the day. 

Figure 19 shov/s values of static stability '■(32-46), 7(24-32), '(16-24) 
as a function of time of day for tlic tower in the clearing. The distribu- 
tion of temperature and stability at those lieiglits are not as strongly 
distorted by the forest canopy as immediately above the canopy, hence, the 
values (and tlie ordinate scale of the diagram) are quite different than in 
the previous diagram for the forest tower stability distribution. The 
stability values t ( 32-46), '(24-32) undergo less diurnal variation than 
a (16-24) and are generally within the range -2 to +5x10 ^m ^ throughout 
the entire 24 hour period. They generally show conditions nearest neutral 
during the day. 

Figure (20) sliows a distribution of values of static stability between 
levels (4-8m) , (2-4m) and (l-2m) along the tower in the clearing. The 
layer next to the surface (l-2m) shovv’s tlie largest diurnal variation of 




STATIC STABILITY (IQ-^m-') 


FOREST. DOTS ARE FOR THE LAYER 36m TO A6m; 
CIRCLES ARE FOR THE LAYER 32m TO 36m; CROSSES 
ARE FOR THE LAYER 25m TO 30m. 




LOCAL TIME 


Figure 18. Diurnal Variation of Static Stability for Various 
Layers Above the Forest Cmopy for January 
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CLEARING. CIRCLES ARE FOR THE LAYER 32m TO A6m; 
DOTS ARE FOR THE LAYER 2Am TO 32m; CROSSES ARE FOR 
THE LAYER 16m TO 2Am. 



Figure 19. Diurnal Variation of Static Stability for Various 

layers in the Clearing for January i 
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FIGURE 20. DIURNAL VARIATION OF STATIC STABILITY AT LOW LEVELS IN THE 
CLEARING FOR JANUARY. 


2 A 


1 J 


-1 


! 

1 

1 

4 

- -4 - - - 

w. 


o o(4-8 METERS) 
• a (2-4 METERS) 
A a(l-2 l^ETERS) 


“•* I‘ 


i 

A 



' ^ A • 

O 0 c * o o 0 t o c 

0 0 ^ 0^0 1 
^ ^ • X 


^ ^QftOOOOC^ 


X - A 


-3-1 


I II I I I 1 II I I II ' i • I I "I I I I ' I I 
o2 4' ^ g /o /J /*/ /(, /g XV 


L <x ca/ 7~i 


93 


stability with a regime strongly coupled to the diurnal variation in 
surface temperature. In the layers (4-8m) and (2-4m) centered at 6m and 
3m respectively, the stability is more uniform over the day although 
j(2-A) exhibits a diurnal cycle similar to o(l-2) but smaller in amplitude. 

Fig. 21 shows a time-height section for two week averaged static 
stability in the forest for January. The section covers a typical 
diurnal variation from the ground to the top of the tower and shows the 
occurrence of unstable, near-neutral, stable and inversion conditions. 
Specifically, the criteria used in this plot is 





(i/e) (dO/oz) 
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X 
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m 
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wtiere = potential temperature 

These criteria correspond rougiily to ttie following 
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As can be seen, near neutral conditions during this period were quite 
rare. The layer of air next to the forest floor shows strong inversion or 
isothermal conditions during the daylight hours from 10:00 to 19:00, then 
lapse conditions set in after sunset and develop large values during the 
night and early morning. The trunk space region is dominated by inversion 
and/or stable conditions with unstable conditions developing below tree-top 
level just before noon and continuing until mid-afternoon. 

After sunrise, a warm layer centered on about the 30m level just above 
the tree tops develops and persists till the mid-afternoon gradually sub- 
siding by sunset. The lower portion ('-25- 30m) of this layer maintains inver- 
sion conditions while the upper portion (—lO-lbm) exhibits lapse conditions 
which become unstable around 10:00. The above canopy jet tends to have its 
maximum speed in the portion of this layer showing significant lapse. During 
the night-time hours, the temperature at 30m shows a gradual cooling which 
establishes a thin unstable layer beneath and an inversion immediately above 
with a secondary unstable layer above that. At higher levels of the tower, 
stable or inversion conditions persia; through virtually the entire day. 
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B. Wind Direction During January and June 

In Part 3 of this Section, It was shown that the wind profile parameters 
- particularly the roughness parameter - varied considerably between 
different periods of the year. There are several factors which might con- 
tribute to these seasonal differences besides the prevailing wind speed as 
discussed in Part 3. One of the most prominent differences in the forest 
environment between January and June is the difference in the direction 
of the prevailing wind. June is in the middle of the southwest monsoonal 
flow while January is characterized by flow from the north or northeast. 

The upwind terrain characteristics during these two months is different 
with several cleared areas about 800-1000 meters northeast of the tower 
in the forest but generally unbroken forest area to the south and southwest. 
Also, the terrain slopes upward from the south and southwest to the towers 
with an inclination of about 1:40 or greater (ASRCT, 1969). The winds 
during June are generally stronger than during January but the direction 
is more variable in January than June. It was noted in Section 11 that the 
wind direction data was logged onto the "A-tapes" and that these tapes were 
not processed during this project. Nevertheless, limited wind direction 
information was available to us on computer printouts made from several 
A-tapes during the TREND field experiment. These are summarized in Table 
20 which shows the diurnal variation of wind direction for January 4, 

1970 and June 23, 1970. The wind direction is in degrees from north. 
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Table 20. 

Wind Direction (Degrees from North) 
and June 23, 1970 

for January 

4, 1970 

• of Day 

January 4, 

1970 

June 23, 

1970 

0000 

47 


262 


0100 

49 


265 


0200 

63 


258 


0300 

102 


259 


0400 

344 


263 


0500 

341 


261 


0600 

346 


258 


0700 

351 


259 


0800 

13 


2 60 


0900 

58 


249 


1000 

53 


247 


1100 

45 


256 


1200 

59 


252 


1300 

30 


253 


1400 

36 


241 


1500 

47 


2 45 


1600 

69 


234 


1700 

83 


254 


1800 

87 


262 


1900 

65 


259 


2000 

42 


2 54 


2100 

76 


270 


2200 

88 


274 


2300 

84 


257 
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C. Spectral Analysis of Wind Speed Data 

A limited amount of spectral analysis was done on wind speed data from 
tlie towers. The ultimate goal is to obtain spectral analyses (of kinetic 
energy and wind-temperature co-spectra) covering a large range of frequencies 
by using half-minute data for the highest frequencies and half-hourly 
averaged data for intermediate frequencies. These studies were started 
relatively late in the contract period and, thus, only kinetic energy spectra 
involving the lower frequencies are presently available. 

The first step in tliese analyses was to develop an appropriate computer 
program for application to the TREND data. The Graduate Program in 
Meteorology has a spectral analysis routine using the Fast Fourier Trans- 
form (FFT) that has been adapted to a number of meteorological studies. 

Results using the FFT were compared with results of direct spectral estimates. 

The time series chosen for analysis with the FFT were for 10 days in 
January and about 10 days in June using data from several levels along 
both towers. The FFT divides the time series into a number of equal and 
overlapping segments, obtains the spectral estimate of eacli segment and 
then averages. The length of the time series was 180 data points (that 
is, 180 30-mlnute averaged values). Since the method requires t'nat the 
lengtli of record sliould be a power of 2, tlie record is supplemented bv zeroes 
up to 256 points whicli is the first power of 2 beyond 180. The number of 
lags was 127 (128 is tlie largest power of 2 which is less than 180). Using 
half-hourly averages as data, 180 values corresponds to 90 hours so that 
tlie largest frequency is about 3 days with this sample wliile the Nyquist 
frequency is (l/2At)=l cycle/hour. 
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The kinetic energy spectral results from the FFT method are shown in 
Figure 22 for the tower in the forest for the period June 20-30, 1970. A 
"red" spectrum was computed as a first order linear Markov process based 
on one lag value of the autocorrelation series. This approach seems 
justified since the two lag value of the series approximates the one lag 
value squared [Gilman, Fuglister and Mitchell (1963)]. The autocorrelation 
values for levels 16m and 46m are plotted in Figure 23. In Figure 22 
two "red" spectra are plotted to illustrate the difference in the manner by 
which the "red" spectra is defined. The dashed line was obtained using two 
lags in the autocorrelation function while the dot-dashed line was obtained 
using one lag. The chi-square test can be used to test the significance 
of each peak. 

For the same record and using the same number of lags (127) , direct 
spectral estimates were obtained and compared with spectra obtained from 
the FFT method. The results are shown in Figure 24 and exhibit features 
comparable to those of Figure 22. That portion of the spectrum above the 
frequency of 20 cycles/hour seems to follow a (-2) power law [and not the 
(-5/3) power law based on the Kolmogorov hypothesis]. 

In order to reduce computation and plotting time and increase tlie 
"stability" of the spectra, the number of lags was reduced to 50. The 
results are shown in Figures 25a, b and are similar in structure to Figure 24. 
On Lne basis of these results, it seems reasonable to deal with the half- 
hourly data by using direct Fourier transforms with 50 lags. n\e results 
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Figure 22. Noriiali/cd Spectral Density of Kinetic Cnergy 

at 4f)m Along the To\;cr in tlie Forest, for June 20- JO, 1970 
Using the Fast Fovirier Tre.nsfons Method t'ii.h 127 i,,ags. 
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FIGURE 25b. NORMALIZED SPECTRAL DENSITY OF KINETIC ENERGY AT 16 METERS 
ALONG THE TOWER IN THE FOREST FOR JUNE 20-30, 1970 USING DIRECT FOURIER 
TRANSFORMS WITH 50 LAGS. 
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comprising Figure 25a were also plotted using the product of spectral density 
and frequency along the ordinate as is sometimes the custom. This display 
is shown in Figure 26. 

Similar computations were done for data taken along the tower in the 
forest during January at levels 46m, 32m, 16m, 4m. The autocorrelation 
functions are plotted in Figure 27. It is readily evident that the levels 
within the canopy (4m, 16m) have different autocorrelation characteristics 
than those levels above the canopy (32m, 46m). The values drop off much 
faster and the computation of the red spectra using one lag in the auto- 
correlation function seems justified. Normalized spectral density for the 
16m level is shown in Figure 28a. In the frequency range above 0.3 cycles/ 
hour, the (-5/3) power law is indicated. Normalized spectral density 
for the 46m level is plotted in Figure 28b. .-t seems that the red spectra 
as computed here with one lag is inappropriate and that in this case it 
would bo better to average over a few lags and then compute the red spectra. 

Before preparing additional computations of spectral estimates, 
the variability of the time series was studied to test one of the basic 
hypotheses - that the data belong to the same population. Using the record 
for June, the data were divided into 5 samples of two days each. The mean 
and standard deviation of each sample were computed as well as the 10 day 
mean. Confidence limits for the mean of tlie population of each sample 
were computed using the "t"-test: 

X - t_Q3 [S/(n-l)^/^] < u < X + t_^3 [s/(n-l)^^^] 




SPECTRAL DENSITY X FREQUENCY 






NORMALIZED SPECTRAL DENSITY 






NORMALIZED SPECTRAL DENSITY 




where x is the sample mean, S is the sample standard deviation, p is the 
population mean, n-1 is the number of degrees of freedom. The following 
results were obtained: For 95 degrees of freedom, t=1.99 and the confidence 
intervals for 95% level of significance appropriate to the 46m level of 
the tower in the forest are shown in Table 21. The ten day average for 
level 46m above the forest was u=4.28 m-sec ^ which falls within the confi- 
dence limits of the first four samples but not the fifth. The confidence 
interval for 98% level of significance for sample 5 is 2.448 < p < 3.072 and 
the confidence interval for the 99.9% level of significance is 2.311 
< p < 3.209 both of which exclude the 10 day mean. 

Similar tests at other levels in both the forest and clearing during 
both June and January are summarized in Tables 22-35 . On the basis of 
these results, spectral estimates for the data taken along the tower in 
the forest during June for a shortened period of 8 instead of 10 days 
were computed. The autocorrelations are not much different, however, and 
are shown in Figure 29. 

The autocorrelations and spectral estimates for both periods for 
the tower in the clearing were computed and are presented in Figures 30-34. 







TABLE 21. Confidence Intervals for the 95% Level of Significance: 

Wind Speed Data at 46m Along the Tower in the Forest 
for June 20-30. 1970 


Level 46m - Forest - June. M = 4.28 m-sec 


Sample Number 


Confidence Interval 


Criterion Satisfied (x) 


Sample 1 


4.121 < u < 4.619 


Sample 2 


3.968 < p < 5.022 


Sample 3 


3.895 < y < 4.405 


Sample 4 


4.210 < y < 4.770 


Sample 5 


3.243 < y < 3.957 


TABLE 22. Same as the Previous Table Except Appropriate to the 40m Level 


Level 40m - Forest - June. y = 3.72 m-sec 


Sample Number 


Confidence Interval 


Criterion Satisfied (x) 


Sample 1 


3.564 < y < 3.996 


Sample 2 


3.984 < y < 4.356 


Sample 3 


3.377 < y < 3.823 


Sample 4 


3.641 < y < 4.119 


Sample 5 


2.866 < y < 3.474 


TA.BLE 23* Same as the Previous Table Except Appropriate to the 36m Level 


Level 36m - Forest 
Sample Number 
Sample 1 
Sample 2 
Sample 3 
S amp le 4 
Sample 5 


June, u = 3.24 m-sec 
Confidence Interval 
3.092 < u < 3.468 
3.461 < u < 3.779 
2.954 < li < 3.326 
3.174 < p < 3.586 
2.499 < u < 3.021 


Criterion Satisfied (x) 

X 

X 

X 


TABLE 24. Same as the Previous Table Except Appropriate to the 4m Level 
Level 4m - Forest - June, p = 0.43 m-sec ^ 


Sample Number 
Sample 1 
Sample 2 
Sample 3 
Sample 4 
Sample 5 


Confidence Interval 
0.435 < p < 0.505 
0.465 < p < 0.535 
0.411 < p < 0.489 
0.345 < p < 0.415 
0.307 < p < 0.373 


Criterion Satisfied (x) 


X 


X 
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TABLE 25. Confidence Intervals for the 95% Level of Significance: 

Wind Speed Data at 46m Along the Tower in the Forest for 
January 4-14, 1970 


Level 46m - Forest - 

Sample Number 
Sample 1 
Sample 2 
Sample 3 
Sample 4 
Sample 5 


January, u = 2.28 m-sec 

Confidence Interval 
2.229 < p < 2.491 
2.744 < p < 2.956 

1.802 < p < 2.018 

2.352 < p < 2.568 
1.734 < p < 1.946 


Criterion Satisfied (x) 

X 


TABLE 26 . Same as Previous Table Except Appropriate to the 32m Level 


Level 32m - Forest - January, p = 0.55 m-sec 


-1 


Sample Number 
Sample 1 
Sample 2 
Sample 3 
Sample 4 
Sample 5 


Confidence Interval 
0.435 < p < 0.565 
0.449 < p < 0.592 
0.307 < p < 0.373 
0.519 < p < 0.662 
0.745 < p < 0.875 


Criterion Satisfied (x) 
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TABLE 27. Same as the Previous Table Except Appropriate to the 4m Level 
Level 4m - Forest - January, u = 0.27 m-sec ^ 


Sample Number Confidence Interval Criterion Satisfied 


Sample 

1 


U 

= 

.2 


Sample 

2 

0.305 

< 

P 

< 0.375 


Sample 

3 

0.256 

< 

P 

< 0.305 

X 

Sample 

4 

0.269 

< 

P 

< 0.331 

X 

Sample 

5 

0.240 

< 

P 

< 0.280 

X 


TABLE 28. Confidence Intervals for the 95% Level of Significance: 

Wind Speed Data at 46m Along the Tower in the Clearing 
for June 20-30, 1970 


Level 46m - Clearing - June, p = 5.9 m-sec 


Sample Number 
Sample 1 
Sample 2 
Sample 3 
Sample 4 
Sample 5 


Confidence Interval 
5.623 < y < 6.267 
6.333 < p < 6.987 
5.294 < p < 6.086 
5.788 < p < 6.572 
4.531 < P < 5.549 


Criterion Satisfied 

X 


X 


X 


(x) 


(x) 
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TABLE 29. Same as the Previous Table Except Appropriate to the 32m Level 
Level 32m - Clearing - June, y = 5.13 m-sec ^ 


Sample Number 
Sample 1 
Sample 2 
Sample 3 
Sample 4 
Sample 5 


Confidence Interval Criterion Satisfied (x) 

4.810 < y < 5.390 x 

5.574 < y < 6.126 

4.655 < y < 5.325 x 

4.983 < y < 5.637 x 

3.949 < y < 4.831 


TABLE 30 . Same as the Previous Table Except Appropriate to the 16m Level 
Level 16m - Clearing - June, y = 3,76 m-sec ^ 

Sample Number Confidence Interva l Criterion Satisfied (x) 


Sample 1 

y = 3.468 


Sample 2 

4.190 < y < 4.590 


Sample 3 

3.475 < y < 3.965 

X 

Sample 4 

3.595 < y < 4.065 

X 

Sample 5 

2.845 < y ^ 3.555 
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TABLE 31 . Same as the Previous Table Except Appropriate to the 2m Level 


Level 2m - Clearing - June, u = 2.35 m-sec 


Sample Number 
Sample 1 
Sample 2 
Sample 3 
Sample 4 
Sample 5 


Confidence Interval Criterion Satisfied (x) 

2.195 < p < 2.465 x 

2.691 < y < 2.949 

2.253 < u < 2.567 

2.191 < y < 2.509 x 

1.591 < y < 2.109 


TABLE 32. Confidence Intervals for the 95% Level of Significance: 

Wind Speed Data at 46m Along the Tower in the Clearing 
for January 4-14, 1970 


Level 46m - Clearing - January, y = 2.95 m-sec 


Sample Number 
Sample 1 
Sample 2 
Sample 3 
Sample 4 
Sample 5 


Confidence Interval Criterion Satisfied (x) 


3.015 

< 

y 

< 

3.305 


3.650 

< 

p 

< 

3.870 


2.731 

< 

y 

< 

3.010 

X 

2.683 

< 

p 

< 

2.957 

X 

2.028 

< 

p 

< 

2.273 
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TABLE 35. Same as the Previous Table Except Appropriate to the 2m Level 
Level 2m - Clearing - January, u = 0.83 m-sec 


Sample Number 
Sample 1 
Sample 2 
Sample 3 
Sample 4 
Sample 5 


Confidence Interval 
0.792 < u < 1.008 
0.974 < y < 1.206 
0.584 < y < 0.816 
0.676 < y < 0.904 
0.554 < y < 0.726 


Criterion Satisfied (x) 
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Figure 32b. Nonr.alized spectral density of kinetic energy at 
32 meters along tlic tower In the clearing for dune 20-30 
1970 using direct Fourier transforms witli 50 lags. 
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RED SPECTRUM (one lag) 



Figure 34b. Normalized spectral density of kinetic energy at 

32 meters along the tower in the clearing for January 4-14, 
1970 using direct Fourier transforms with 50 lags, 
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SECTION V 
ENERGY FLUXES 


1. Radiatlonal EnerRy 
A. Data Available 

As noted in Section II of this report, radiation and dew point data 
from the TREND field experiment was logged onto A-tapes and none of these 
data were preprocessed into a readily usable form. Some limited amounts 
of radiation and dew point temperature data were available in the form of 
computer printouts of A-tape data made during the TREND experiment. Thus, 
extensive analyses of the radiation (and latent heat) budgets of the Thai 
forest area are not possible at this time. Representative samples of 
radiational data are presented in this section using the printout informa- 
tion at our disposal. 

Table 36 gives a summary of the periods for which measurements of 
radiation fluxes, wind direction and dew point temperature could be taken 
from computer printout sheets. As noted in this table, (and in Section II), 
much of the recorded dew point temperature information is faulty due to 


malfunctioning sensors. 
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TABLE 36. 


Log of Periods for Which Radiation, Wind Direction and Dew 
Point Temperature Measurements are Available from Computer 


Printouts Made During the TREND Field Experiment 


Radiation 


Wind Direction, Dew Point Temperature 

February 6-23, 25-28 

1969 

*February 6-15, 18-23, 26-28 

March 1-13, 18-25 


*March 1-4, 23-25, 27-31 

April 1-13, 20-26 


*April 1-5 

September 29-30 



October 1-10, 11-12 


*0ctober 28-31 



*November 1-6 

December 22-24 


*December 22-24 

January 7-19 

1970 

^January 3-21 

February 5-24 


*February 5-11 

March 3-8, 14-29 


*March 4-13 

April 5-14, 20-30 


*April 25-30 

May 1-18, 20-31 


*May 4-5, 27-31 

June 1-4, 15-22 


June 1-13, 22-30 

August 14-22, 27-31 


August 1-4, 25 

September 4-18, 21-29 


*September 4 


* Faulty Dew Point Data. 




136 


B. Albedo Determinations 

With a knowledge of incoming solar radiation one can estimate the 
outgoing stream by multiplying the incoming flux by the overall albedo, 
if known. Using measurements of both and Rg*i one can estimate the 
albedo from the ratio [R^t/R^+l. With a representative value for the albedo 
for the forest, we reduce by one the number of radiative fluxes that must 
be measured for an energy balance survey. Table 37 shows estimates of 
albedo made from radiation measurements during the daytime in each of four 
different months of the year. The overall average value of albedo for 
these determinations is about 14% but varies from about 11% to 17% in the 
sample tested. There appear to be some consistent differences with season 
in the sample. Since the forest top is not a rigid surface, one expects 
variations in albedo related to the environmental conditions in the 
forest. 


C. Net Radiation 

Table 38 and Figures 35-38 represent tlie information on the net radia- 
tion in a number of characteristic days during spring, summer and winter. 
December 23, 1969 and April 21, 1969 are clear days and seem to be perfect 
for energy balance study; however, no humidity data were available for these 
days. The required corrections to the incoming and outgoing radiations 
were made for the August 11-13, 1970 data, since the Funk black body 
temperatures necessary to make these corrections were on the printouts. 


table 37. 

DATE 

3/21/70 

6/15/70 

9llin^ 

\inilb9 


Albedo Estimates for a Tropical Dry Evergreen Forest for 
Selected Seasonal Periods 


TBfE R + , , -7 -1 albedo 

(ca?-cm -min ) (cal-cm -min ) (%) 


0930 

0.7420 

0.0^196 

13.4 

1030 

0.8663 

0.1065 

12.3 

1330 

1.0370 

0.1153 

11.1 

1430 

0.8853 

0.1046 

11.8 

1530 

0.3959 

0.0497 

12.6 









1630 

0.0927 

0.0180 

19.4 

0930 

0.2734 

0.0346 

12.7 

1030 

0.5203 

0.0654 

12.6 

1330 

0.5088 

0.0614 

12.1 

1430 

0.6902 

0.0847 

12.0 

1630 

0.3297 

0.0456 

13.3 




12.6 





1730 

0.0806 

0.0144 

17.8 

0930 

0.5594 

0.0863 

15.4 

1100 

0.8424 

0.1150 

13.7 

1400 

1.2116 

0.1408 

11.6 

1500 

0.7195 

0.0854 

11.9 

1600 

0.1164 

0.0195 

16.7 




13.9 


0930 

0.6975 

0.1215 

17.4 

1030 

0.9227 

0.1481 

16.0 

1330 

0.9982 

0.1539 

15.4 

1430 

0.8092 

0.1295 

16.0 


•Average 16.2 


Overall Average 14.1% 
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The numbers In parenthesis represent 


i 


the difference between corrected and 
within the limits of possible error. 


the corrected values. It seems that 
uncorrected net radiation values is 


TABLE 38. Net Radiation (Langleys/day) for Several Days During the 
TREND Experiment 


n 

DATE 

NET RADIATION 

% 

FF/FT 

Hours of j 

Sunshine 

Forest top 
(Ly/d) 

Forest floor 
(Ly/d) 

April 21, 1969 

441 

43 

9.75 

10.5 

August 11, 1970 

288 (288.6) 

16 (16.7) 

5.5 

> 8 

August 12, 1970 

306 (306.9) 

21 (20.1) 

00 

1 

> 8 

August 13, 1970 

287 (287) 

15 (14.5) 

5.2 

> 8 

August 30, 1970 

317 

29 

9.1 

CNl 

O 

December 23, 1969 

353 

-18 



8.9 






DIURNAL M^IATION OF NET RADIATION 2iaprili969 




NET RADIATION ii august isto 



Figure 36. Diurnal varialion of net radiation above tiie forest 
and at the lorest floor for August 11. 1970. 



Figure 37. Diurnal variation of net radiation above the forest 
and at the forest floor for A\igust 30, 1970. 
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Some attempts were made to parameterize the net radiation flux in 
terms of various temperatures measured simultaneously along the towers. 

For energy balance studies, direct measurements of radiation are, of course, 
much more desirable than any parameterization but since the record of 
radiational data is quite limited, it is of interest to search for empirical 
estimators of such quantities. It is also true that a parameterization of 
radiational fluxes should incorporate information on cloud cover, humidity, 
etc., but once again, the motivation was to test parameterizations involving 
only temperature since there is far more temperature and wind speed in- 
formation for these studies than any other kind. 

The initial test was to search for a simple correlation between the 
net infra-red radiation flux and the temperature in the forest canopy, 

T^(“K) [actually, the estimator used was a=Stephan-Boltzmann constant]. 

Radiation and temperature data were selected during nighttime hours between 
2230 and 0530 for the dates June 19, 20, 21, 22, 1970. For the resulting 
sample of 45 measurements we found the regression relation as follows: 




= 0.1777(aT^) - 0.1047 
c 


where the radiation is in cal-cm -min , is the absolute temperature 

at the 30 meter level of the tower and is taken to represent the canopy 

temperature. The correlation coefficient was only 41.9% and the standard 

-2 -1 

error of estimate was 0.00178 cal-cm -min . Since this correlation was 



y 
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rather weak, the data were submitted to a stepwise screening regression 
algorithm with several (admittedly interdependent) temperature predictors 
to discover which were the most significant predictors and how great a 
correlation could be obtained using temperature predictors alone. The 
predictors used were as follows: 

4 

P^=(aT^) ; T^= canopy temperature [temperature at the 30 meter level] 

4 

P^=(aT^) ; T^= temperature at the forest floor 
4 

P =(oT ) ; T = average temperature for the layer 0 m to 46 m height 

a ^ cL 

4 

Pj.= (aT^) ; T^= average temperature for the top layer 30m to 46m height 



During the nighttime hours for clear sky conditions, the net infra- 
red radiation is generally negative; that is, more energy is radiated 
upward by the surface and/or forest than is received from the atmosphere 
above. The first variable selected by the stepwise regression scheme was, 
consequently, the forest floor temperature, T^ , which turned out to be a 
more powerful predictor than the canopy temperature. Of the four temperature 
predictors mentioned above, the above-canopy average temperature, T^ , should 
be most representative of the radiating characteristics of the atmosphere 
above the forest and, hence, of the Incoming infra-red radiation. Indeed, 
the second predictor selected by the stepwise regression sclieme was T^ . The 

third and fourth predictors were, respectively, iT and T . The regression 

c 

equations corresponding to the four steps of the stepwise regression scheme 


are given in Table 39 . Our interpretation of these results is that 


i;: VARIOUS layurs of ti:e forest 
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temperature quantities alone are not sufficiently powerful as predictors 
and that one should select only and and test other non-temperature 
predictors (relative humidity, cloud cover, etc.) for their ability to 
reduce variance. Also, the constant terms in each of the equations above 
are rather large which indicates that the time mean value of the (radiating 
temperature) predictors are rather significant in determining the net 
radiation. In fact, this has a great deal to do with the resulting algebraic 
signs which precede the various predictor terms in Table 39 . While the 
regression equations could be linearized and the means extracted, one cannot 
be too encouraged at the prospects of parameterizing radiational fluxes in 
terms of temperature alone especially when one recalls the large variability 
of cloudiness and humidity during the various monsoon periods. 

2. Sensible and Latent Heat Fluxes 

A. Outline of Possible Estimating Procedures 

One of the most interesting aspects of the TREND experiment was the 
potential for estimating all components of the complete energy and moisture 
balance of a tropical evergreen forest for many periods during the complex 
monsoon cycle. As the reader is by now aware, it was not feasible to 
completely achieve this goal during this phase of the tropical environmental 
study due largely to the rawness of the original data and certain sensor 
malfunctions. VJith the large amount of temperature and wind speed data, 
it was possible to make some estimates of sensible heat flux and limited 
latent heat flux estimates. Before presenting these results, it is of 
interest to review briefly several alternative procedures for making such 


estimates . 
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Gradient Method 


One of the better known methods for indirect flux determination 
utilizes the gradient of temperature, specific humidity and velocity as 
expressed in the following relations : 


H = - CpPK^(dT/dz) (1) 

E = - pK (dq/dz) (2) 

w 

T = pK^(du/dz) ( 3 ) 

where 

H - sensible heat flux (ly/min) 

3 

E - rate of evaporation (gr/cm sec) 

T - shearing stress 

3 

p - air density (g/cm ) 

2 

- eddy diffusivity for horizontal momentum (cm /sec) 

Cp - specific heat of air (cal/gdeg) 

2 

- eddy diffusivity for heat (cm /sec) 

Equation (3) is introduced here to point out the assumption of similarity 
between the vertical flux equations of momentum, heat and moisture. 

In order to bring the previous relations into operative status it is 
the practice to make some assumptions and modifications about the relation 
of the turbulent transfer coefficients. According to Pasquill (1949) they 
are Identical only under stable and neutral conditions. Swinbank (1955) 
has shown that K^/K^ Increases with increasing instability. Taylor (1960), 
after analyzing Swinbank' s data and the data of Rider (195y ) comes to the 
conclusion that both factors are identical. Uncertainty remains only under 
very unstable atmospheric conditions. Investigations by Crawford (1965) 
favor the identity of both factors. To employ this method, one must make 


direct use of estimated values of di f fusivities . 
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Bowen Ratio Method (Bowen, 19 26) 

One of the methods which avoids direct use of the eddy diffusivity 
factors but uses only their quotients is the Bowen ratio method. The Bowen 
ratio is defined as: 

T - T 

B = H/LE = Y (K. /K ) ( ) 

li w e -a 
s a 

where 

L - latent heat of cv^aporation 
T ,T - surface and air temperature 

S 3l 

e^,e^ - vapor pressure at surface and in air 
, - psychrometer constant (0.485 mmHg • °C 

Assuming identity of transfer coefficients VL , K and using the heat budget 

11 w 

equation adjusted to the air/forest interface, one may write, 


R + S + EL + H = 0 


where 


R - the net radiation at top of forest 

S - heat gained or lost by vegetation soil and photosynthesis 
EL,H - as before 

If B = 1/6, the energy balance equation may be solved to yield 
Q = - (R + S)/(l + B) 


EL = - (R+S)B/(1 + B) 

This method has been used by many investigators and has been shown to 
give reliable measurements in many agricultural situations (Tanner, 1968; 
Albrecht, 1937; Sverdrup, 1951 ), 
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Profile Method 

Another classical approach to estimating fluxes of heat and moisture 
is assumption of adiabatic conditions and integration of the following: 


where 


K = K = ku*z 
w ra 


k - von Karman constant 
u* - friction velocity defined as 
u* = (t/p)^^^ 


Combining (3), (4) and (5) yields 


* = ^^^2 - ^ 1 ^ 

^ ln[(z„-d+z )/(z -d+z )] 
z o 1 o 


(4) 


(5) 


( 6 ) 


where 

Uj^, U 2 are wind speeds measured at two heights above the canopy z^ and 
Z 2 , and d is the zero plane displacement (explained later). Substituting 
(6) into (4) and then into (2) yields 


E = 


-pk (U 2 - (q^ - 

z„-d+z _ 
1 o 


which is known as the Thornwalte and Holzman equation. If (6) and (4) are 
substituted into equation (’) instead of (2), one obtains 


-CpDK^(u, - u^)(T2 - T^) 

z,-d+z _ 

1 o 



H 
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Eddy Correlation Method 

Vertical transport of heat and water vapor in the atmosphere is mostly 
turbulent. If we define q' as the instantaneous deviation of q from its 
mean value at the point, the product q'w'(w' - the vertical component of 
the wind speed) is the deviation of the instantaneous water vapor flow 
from its mean value. The evaporation is found by averaging over time the 
quantity 

E = cw'q’ 

Similarly for sensible heat flux: 

H = C tw'T* 

P 

Where the overbar represents the averaging operation. 

This best method of estimation is difficult to apply to the TREND 
data since there are no direct observations of vertical velocity and the 
only values of specific iiumidity available are values related to tliirty 
minute averaged values of dewpoint temperature. 

B. Results 

From Table 3b it is clear that very little dew point temperature data 
is available for computational purposes. [A large amount of the dewpoint 
temperature data appearing in the computer printouts of A-tape data are 

_3 

values of order 10 °C ! ] Values of Bowen's ratio were computed for four 
days in June using data at the 32m and 36m levels along the tower in the 
forest which was not obviously and absurdly faulty by virtue of its magnitude. 
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Computations were made for each half-hour period and then averaged over 
six-hour periods of each day. Resulting values are shown in Table 40. 

TABLE 40. Bowen's Ratio for Four Days in June, 1970 


Bowen Ratio Time of Day 


Date 


0:00-6:00 

6:00-12:00 

12:00-13:00 

13:00-24:1 

J une 2 3 , 

70 

.3729-03 

-.1468-01 

.4620-02 

.9214-03 

June 24, 

70 

-.2253-02 

-.7589-02 

-.5274-02 

-.5352-03 

June 25, 

70 

.6290-03 

-.1788-01 

-.4522-02 

.1124-02 

June 26, 

70 

.6101-03 

.8557-03 

.1583-03 

-.7583-03 


The values shown in Table 40 are extremely small and fluctuate in sign. 
According to the stability statistics given in Section IV, the layers above 
the forest canopy are frequently near neutral conditions so tiiat a Bowen 
ratio would be the ratio of two small and variable energy fluxes. It is 
not clear that the ratios appearing in Table 40 are very revealing. To 
further document the situation on these four days, however, six-hour 
totals of rainfall at the top and floor of the forest were computed. The 
rainfall totals are shown in Table 41 where it may be seen that rain indeed 
fell during each of the 16 six-hour periods. Also, the fraction of the 
daily total rainfall which reached tlie floor of the forest varied from a 
high of 93% on June 26, 1970 to a low of about 49% on June 23, 1970. The 
rain gauge at the forest floor actually collected more water than fell as 
rain for 3 of the six-hour periods shown in Table 41. 
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Having performed determinations of roughness parameter and datum 

displacement level for a large sample of wind profiles during January, 

June and September, the sensible heat flux [and, perhaps, the latent heat 

flux] may be calculated using the Thornthwaite-Holzman equations given 

in the previous subsection. Such calculations were carried out for the 

periods in January, June and September noted in Table 17 and will be 

presented shortly. Before presenting these results, attention is drawn to 

the fact that there was a wide variability in the values of z and d 

o 

calculated using Steam's method for the three periods, liliile the Thorn- 
thwaite-Holzman equation is linear in the wind and temperature information, 
it is highly non-linear in z^ and d. In computing, say, sensible heat 
flux from tills equation, one faces a choice of using some average, repre- 
sentative values of z^ and d in the equation or using "instantaneous" 
values of these parameters determined from the actual profile from which 
the heat flux is computed. If the variability of these "site parameters" 
was sufficiently small, one would feel reasonably comfortable in using 
average values in heat flux computations. Table 17 of Section IV gives 
evidence that such is not the case here. This issue is quite relevant to 
micrometeorological studies since some methods for estimating z^ and d 
yield only one pair of estimated values and the investigator has no choice 
but to use this pair in all heat flux computations. 

One can get an estimate of how bad the situation might be by performing 
an error analysis of the lieal flux equation of Thornthwaite and Holzman. 
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TABLE 

41. Rainfall 
for Four 

at Forest Top (FT) 
Days in June, 1970 

and Forest 

Floor (FF) 





Rain- 





Daily 


Date 

fall (mm) 


Time 

of day 


Total 





0-6 

6 -12 

12-18 

18-24 


June 

23, 

70 

FT 

6.6 

5.7 

6.0 

6.9 

25.2 




FF 

2.4 

3.0 

3.6 

3.3 

12.3 

June 

24, 

70 

FT 

4.8 

5.7 

3.3 

6.9 

20.7 




FF 

2.7 

3.0 

2.7 

3.3 

11.7 

June 

25, 

70 

FT 

7.5 

2.1 

3.0 

4.2 

16.8 




FF 

2.7 

2.7 

3.0 

2.7 

11.1 

June 

26, 

70 

FT 

6.6 

41.1 

6.0 

9.3 

63.0 




FF 

2.3 

41.4 

11.1 

2.7 

58.5 


If one defines AH as the error incurred in calculating the heat flux due 
to errors, Az^ and Ad, in estimating the roughness and datura displacement 
height, then one may write 

AH = (8H/3z ) Az + (3H/3d)Ad 
o o 

The maximum relative error is obtained by dividing both members of this 
equation by H and summing absolute values: 

(AH/H) = |(l/H)(3H/3z ) Az I + I (1/H) (3H/3d) Ad | 

I 0 0 
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The relative derivative terms may be calculated from the heat flux equa- 
tion. Carrying this out one may write 

(AH/H) = C [|(Az^/z^)| + I (Ad/z^) I ] 


where 


C = [2 z^(z 2 -z^)] / Y^ In (Y^/Y^) ] 

Y = z - d + z ; Y = z„ - d + z 
11 o 2 2 o 

Thus, the error effect Involves only z^ , d, and the pair of levels u£;ed in 

the computation. One may estimate the error effect from the variability 

of z^ and d shown in Section IV in the following way: if one uses average 

values z and d for z and d and standard deviations in z and d for 
o o o 

Az^, Ad, then (AH/H) represents a measure of the discrepancy involved in 
calculating H using average rather than instantaneous values of tlie profile 
parameters. That is, the quantity (AH/H) calculated in this manner becomes 
an estimate of the error limits encompassing 2/3 of the values of 
[H(z^, d)-H(z^,d) ]/H(z^,d) . Table 42 gives values of (AH/H) appropriate 
to the statistics of z^, d for the periods in January, June and September 
Summarized earlier. It is seen that the non-linear error effect is actually 
a damping effect. For example, for the sample covering September 9-11, 

1970 using the levels 40m and 46m, the standard deviations of z^ and d 
relative to z^ are 35.4% and 21.3% respectively, yet the coefficient C. = O.J97 
(which is a damping of error since C < 1) and tlie discrepancy (AH/H) = 

22.5%. Thus, even though tt>e errors in and d are additive in Ibis 


TABLE 42. Estimates of Relative Error in Sensible Heat Flux Calculated Using the Thornthwaite- 
Holtzman Eauation with Mean Values of Displacement Height and Roughness Parameter. 




f 

i 


approach, the corresponding relative error in H is less than the error in 
determining alone. It is encouraging that errors in determining 
and d are damped by the Thornthwaite-Holzman equation. Nevertheless, the 
entries in Table 42 suggest that one may incur large discrepancies in 
calculating H using mean values of and d as opposed to instantaneous 
values. The coefficient C is larger in June than September and January 
for all three layers shown and generally increases with decreasing height 
of the layer. Due to larger variability in and d in September compared 
with June, however, the discrepancy in H is larger in September than June. 

The value of C is minimum in January so that even though the variability 
of z^ is largest in this month, the discrepancy in H is smallest. 

To obtain an empirical estimate of the effect of using period mean 

profile parameters, values of sensible heat flux were computed for June 23, 

1970 from the TREND profile data using the Thornthwaite-Holzman formula witli 

period average values of z^ and d, and then using "instantaneous" values of 

z^ and d. The result is shown in Figure 39. The differences are small when 

tne heat flax is small but of increasing importance as the flux becomes l.iri’C. 

Generally speaking, however, the two methods give compatible results. The 

observed differences at each half-hour interval during -lune 23, noirial izeJ t ■ 

the heat flux value calculated with the half-hourly values of z and d. (tli.it 

o 

is, (AH/H) observed] were compared with a theoretical estimate of the Jiscreti- 

ancy calculated as indicated previously but using z , d , Az = z - .z , ind 

o o o o 

Ad = d - d, for each half hour period [that is, (.Ml/H) estim.ite]. Tiie rel i- 
tion between the two is shown in Figure 40. lost of the values lie on i 



straight line indicating that most of the discrepancy between tin- curvi'.s in 
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32-36 meters above the forest canopy computed using the 
Thornthv.'aite-llolzman equation with average and instan- 
taneous values of roughness and displacement height. 


(AH/H) estimated 
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Figure 39 may be attributed to the variation in roughness parameter and dis- 
placement height throughout the period. The slope of the relationship shown 
in Figure 40 is not unity, however, which indicates that the discrepancy be- 
tween the two methods of calculation is actually less than that predicted by 
the theoretical error analysis. 

Finally, Figures 41a, b, 42a, b, 43a, b show values of sensible heat flux 
for two days each of the months of January, June and September, 1970. Compu- 
tations were made using mean values of the profile parameters for the three 
periods for each of three layers above the forest canopy. The characteristic 
temperature structure above the canopy in January, which was discussed in 
conjunction with Figure 10a, yields strong upward sensible 'leat flux from the 
canopy during the day through t.Pe layer 32-36m but variable or downward fluxes 
through the higher layers due to the tendencv for Inverted temperature struc- 
ture in these higher layers. The total transport upward from the canopy was 
larger for January 8, wiien there were nine hours of sunshine falling on the 
canopy than for January 6 when there were only 1 . 1 hours of sunshine on the 
canopy. Tlie sensible heat flux in June tends to be positive through all the 
layers during the day with slight downward flux at night in the layer immediate- 
ly above the forest canopy. Both days were partly cloudy with June 20 ex- 
periencing more than three but less than elglit hours of sunshine at canopv level 
but June 21 had only a half-hour of sunshine on the canopy during the dav. Note 
that the scales of Figures 42a and b are different and that the total transport 
was much smaller on June 21. Botli days had intermittent rainfall which further 
tended to Interrupt the daytime surge of sensible heat flux from the canopv. 

The two days in September for which Figures 43a, b were prepared had little 
sunshine and significant rainfall. Wti 1 e the sensible heat flux tends to he 
upward during tlie daytime hours, tlie values are small and the variabilitv large. 



SENSIBLE HEAT (cal-cm"*-min"') 



Fj^nre Ain. Sensible hcuit flux Lliroufji various layers above 
tlie forest for Jt.r.uary 6,1970. 




SENSIBLE HEAT FLUX (cal- cm"2-min~') 
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LOCAL TIME 


Figure Alb. Sciisiblo beui; flux through various layers above 
the forest for .lanuiry 8, 1970. 


LOCAL TIME 
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o 42a. Sc,nr;iblc liP.at fliLX Lhrou^li vaiioiis layers above 
the for June 20, 1970. 





uuo- 





SENSIBLE HEAT FLUX (cal-cm"2-min"') 



A'ia. SriTiihlo lu'aL flux thfuuj'.h various Jayors ahovo 
Llu’ forust. for fiopl cmhor 1), 1970. 
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SECTION VI - GENERAL SUMMARY 

Tills report has documented the major results of work done at tlie 
University of Maryland over the two year period June, 1972 - August, 1974 
under Contract DAAK02-72-C-0287 to the U. S. Army Engineer Topographic 
Laboratory, Ft. Belvoir, Virginia. The general objective of this work 
has been to extract scientific information from various raw meteorological 
data materials collected during Project TREND (Tropical Environmental 
Data). Project TREND was a large U. S. Army sponsored interdisciplinary 
scientific field program in the Sakaerat Forest of Thailand conducted 
during the years 1967-70. 

Certain ground station meteorological data were collected during the 
period Nov. 1967 to Sept. 1970 and certain other micrometeorological data 
were collected along two towers in the forest area during about the last 
two years of the field experiment. At our disposal were measurements of 
temperature, relative humidity, precipitation, evaporation, hours of 
sunshine, solar radiation, infra-red radiation, profiles of wind, temperature 
and dew point temperature in the forest canopy and in a cleared area within 
the forest area, and profiles of temperature beneath the ground at both 
sites. As the original data tapes logged at the TREND experimental site 
were in very rough form when they were turned over to the University of 
Maryland, a considerable effort was required to develop computer procc’ssing 
techniques which would produce a usable data source for scientific analysis. 
This report documents the condition of the raw data and siunraarizes the 
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current usability and availability of TREND Project micrometeorological 
data. 

Data from the TREND project have been analyzed to yield a coherent 
depiction of the climatological conditions at the Thai forest site. A 
comparison was made of climatological conditions in the Sakaerat Forest 
with conditions at a number of locations in Southeast Asia. Locations chosen 
for comparison generally yield a depiction of the monsoonal cycles on the 
eastern and western coasts and at low, middle and higher latitudes in the 
interior of the Southeast Asian land mass. It is found that the local 
effects of the large scale monsoon system in this area of the world are a 
strong function of orography and windward proximity to the tropical waters 
in which the area is nested and that mesoscale variations in, say, precipi- 
tation may be quite dramatic in Southeast Asia. 

The influence of the tropical forest on the nature of the vertical 
profiles of wind and temperature have been given in this report. Using 
profile data on a frequency of once per thirty minutes, careful estimates 
of forest roughness, datum displacement height and frictional stress have 
been made for three periods of the year: January (during the cool, dry, 
northeast monsoon season); June (in the middle of the warm, moist, south- 
west monsoon season); September (during the month of maximum precipitation, 
also in the southwest monsoon season) . These estimates have been compared 
with similar information for other types of crop canopy. Spectral analyses 


of these data have also been carried out for periods of about 10 days in 
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length during January and June. 

Certain elements of the energy budget of the tropical forest were 
studied. These studies were somewhat limited because of difficulties with 
extracting usable radiation and dew point temperature data from the original 
TREND experiment data tapes and certain instrumental malfunctions during 
the field experiment. Estimates of sensible heat flux were made and are 
presented. Some calculations of latent heat flux and Bowen's ratio were 
made using the dew point temperature measurements from the TREND project. 
Information on incoming and reflected solar radiation fluxes and on the 
albedo of the Sakaerat Forest are presented. Diurnal variations of net 
radiation above the forest canopy and at the floor of the forest are given 
for several days distributed through the year. 

Certain aspects of the research initiated during this project will be 
continued. In particular, attention will be given to the following: 
extending the survey of wind and temperature profile characteristics of the 
tropical forest region to include diabatic conditions; estimation of 
sensibje heat fluxes under near-neutral and diabatic conditions for 
several periods of the year; analysis of other components of the energy 
budget including conductive heat flux into the soil; extension of spectral 
analyses of wind and temperature data to higher frequencies by utilizing 
the once per half-minute profile data; and further attempts to extract 
usable radiation and dew point information from raw A-tapes for various 
periods of the year. Sufficient success in these endeavors would provide a 
comp re liens ive survey of the energy budget of the tropic.il forest region 
which would be of wide interest to environmental scientists. 
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